ASYMPTOTIC EXPANSIONS IN THE CLT IN FREE PROBABILITY 
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Abstract. We prove Edge-worth type expansions for distribution functions of sums 
of free random variables under minimal moment conditions. The proofs are based on 
the analytic definition of free convolution. We apply these results to the expansion of 
densities to derive expansions for the free entropic distance of sums to the Wigner law. 



1. Introduction 

In recent years a number of papers are devoted to limit theorems for the free convo- 
lutions of probability measures. Free convolutions were introduced by D. Voiculescu [40], 
[41]. The key concept is the notion of freeness, which can be interpreted as a kind of in- 
dependence for noncommutative random variables. As in the classical probability where 
the concept of independence gives rise to the classical convolution, the concept of free- 
ness leads to a binary operation on the probability measures, the free convolution. Many 
classical results in the theory of addition of independent random variables have their 
counterparts in Free Probability, such as the Law of Large Numbers, the Central Limit 
Theorem, the Levy-Khintchine formula and others. We refer to Voiculescu, Dykema and 
Nica [42] and Hiai and Petz [25] for an introduction to these topics. Bercovici and Pata 
[11] established the distributional behavior of sums of free identically distributed ran- 
dom variables and described explicitly the correspondence between limit laws for free and 
classical additive convolution. Using subordination functions for the definition of the ad- 
ditive free convolution, Chistyakov and Gotze [22] generalized the results of Bercovici 
and Pata to the case of free non-identically distributed random variables. It was shown 
that the parallelism found by Bercovici and Pata holds in the general case of free non- 
identically distributed random variables (see [13] as well). This approach allowed us to 
obtain estimates of the rate of convergence of distribution functions of free sums. An 
analog of the Berry-Esseen inequality was proved for the semicircle approximation in [22]. 
For related results see [26]. 
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In this paper we obtain an analogue of Edgeworth expansion in the Central Limit Theo- 
rem (CLT for short) for free identically distributed random variables, based on the method 
of subordination functions. In addition we shall give a bound for the remainder term in 
this expansion. In order to deduce this expansion we establish an approximation of dis- 
tribution of normalized sums of free random variables by the free Meixner distributions. 
An interesting feature of our expansions is that they hold without the not moment re- 
lated assumptions used in the classical probability. For bounded free random variables 
we give an asymptotic expansion for the density of normalized sums. With the help of 
this expansion we obtain the rate of convergence in the entropic free CLT. 

The paper is organized as follows. In Section 2 we formulate and discuss the main 
results of the paper. In Section 3 and 4 we formulate auxiliary results. In Section 5 we 
describe a formal expansion in the Free CLT and in Sections 6 and 7 we prove Edgeworth's 
expansion in the CLT for free identically distributed random variables. Since the proofs of 
Theorem 2.1 and Theorem 2.3 (see Section 2) are similar, we give a proof of Theorem 2.3 
in details in Section 5 and an outline of the proof of Theorem 2.1 in Appendix. In Section 
8 we obtain the local CLT for free bounded identically distributed random variables. In 
Section 9 we obtain the rate of convergence in the free CLT in total variation metric and 
in Section 10 we derive the asymptotic expansion for the free entropy of normalized sums. 



Denote by M. the family of all Borel probability measures defined on the real line 1R. 
Define on M the compositions laws denoted * and EH as follows. For v e M., let /i * v 
denote the classical convolution of \x and v. In probabilistic terms, fi*u = C(X+Y), where 
X and Y are independent random variables with \x = C(X) and v = C(Y), respectively. 
Let fi EH v be the free (additive) convolution of /i and v introduced by Voiculescu [40] 
for compactly supported measures. Free convolution was extended by Maassen [30] to 
measures with finite variance and by Bercovici and Voiculescu [9] to the class Ai. Thus, 
\x EH v = C(X + Y), where X and Y are free random variables such that /i = C(X) and 
v = C(Y). There are free analogues of multiplicative convolutions as well; these were first 
studied in Voiculescu [41]. 

Henceforth X, Xi, X 2 , . . . stands for a sequence of identically distributed random vari- 
ables with distribution ji = C(X). Define 



where k — 0, 1, . . . and q > 0. 

The classical CLT says that if Xi,X 2 , . . . are independent and identically distributed 
random variables with a probability distribution /i such that mi = and m 2 = 1, then 
the distribution function F n (x) of 
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X x + X 2 + ■ ■ ■ + X, 
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(2.1) 
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tends to the standard Gaussian law <&(x) as n — > oo uniformly in x. 

A free analogue of this classical result was proved by Voiculescu [39] for bounded free 
random variables and later generalized by Maassen [30] to unbounded random variables. 
Other generalizations can be found in [10], [11], [22], [26]-[28], [35], [45], [46]. When 
the assumption of independence is replaced by the freeness of the noncommutative random 
variables X±,X 2 , . . . , X n , the limit distribution function of (2.1) is the semicircle law w(x), 
i.e., the distribution function with the density p w (x) := 2~a/(4 — x 2 ) + , x e K, where 
a + : = max{a, 0} for a6l. Denote by /i w the probability measure with the distribution 
function w(x). 

Write <f(x) := -J=e~ x2 / 2 and denote by H m (x) := (— l) m e^ 2 / 2 J^e - ^/ 2 the Hermite 
polynomial of degree m. 

Assume that the random variables Xj are independent and have moments of all orders. 
For the distribution function F n (x) of Y n there exists a formal expansion in a power series 
in (see [24], [36]): 

F n (x)^Hx) + <p(x)J2^ (2-2) 

p=i 

where 

and 7 m is the cumulant of order m of random variable X. In the last equality the summa- 
tion on the right-hand side is carried out over all nonnegative integer solutions (ki, . . . , k m ) 
of the equations 

ki + 2k 2 + • • • + pk p = p and s — k± + • • • + k p . (2.3) 

Note that Qi(x) = —m 3 H 2 (x)/6. 

In terms of characteristic functions (2.2) has the form 



oo 00 

2 , 



e to dF n (x) = e-' 2 / 2 + ^§ e ~ t2/2 > (2-4) 

J — oo 

where 



n m/2 

m=l 



oo 

e itx dQ m {x)=P m (t)e- t2 / 2 



Esseen [23] proved that if the random variables Xj are independent, non- lattice dis- 
tributed and f3 3 < oo, then F n (x) admits the following asymptotic expansion 

F n (x) = $(x) - 0=H 2 (x)<p(x) + o(l/ Vn) (2.5) 
which holds uniformly in x. 
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If the random variables Xj are independent and are lattice distributed, that is they 
take values in an arithmetic progression {a + kh; k — 0, ±1, . . . } (h being maximal), and 
f3 3 < oo, then 

F n (x) = <&(*) + ±=*p{x) ( - ^H 2 (x) + hT(p& - ^)) + o(l/v^), (2.6) 

uniformly in x, where T(x) := [x] — x + 1/2. 

If absolute moments (5k of order k > 3 exist, then generalizations of the asymptotic 
expansions (2.5) and (2.6) hold under additional conditions on the characteristic function 
of A" [36]. 

An analytical approach using subordination functions allowed us to give explicit esti- 
mates for the rate of convergence of distribution functions of Y n in the case of free random 
variables. We demonstrated this (see [22]) by proving a semicircle approximation theorem 
(an analogue of the Berry-Esseen inequality [36], p. 111). In this paper we shall establish 
Edgeworth expansion in the semicircle approximation theorem and a complete analogue 
of the Berry-Esseen inequality for identically distributed free random variables. 

We now formulate the main results of the paper. As before we denote by F n (x) the dis- 
tribution function of Y n where Xj are free random variables with the same distribution /i. 
Assume as well that Xj have moments of arbitrary order and m\ = 0, to 2 = 1. We denote 
by /!„ the distribution of Y n . Denote by U m (x) the Chebyshev polynomial of the second 
kind of degree to, i. e., 

U m (x) = U m (cos0) : m = l,2,.... 

It is easy to see U\(x) = 2x, U 2 (x) = Ax 2 — 1, U 3 (x) = Ax(2x 2 — 1). 

It turns out that there exists an analogue of the formal expansion (2.4) for F n (x). To 
formulate it we need the following notation. Define the Cauchy transform of ji G M. by 

GM = [ z E C+, (2.7) 

J R z — x 

where C + denotes the open upper half of the complex plane. The formal expansion has 
the form 

fi^^W + E^, (2-8) 

k=l 

where 

^) = E^ (l//- Z )m ( 2 - 9 ) 

with real coefficients c Ptm which depend on the free cumulants 0:3, . . . , ak+2 and do not 
depend on n. The free cumulants will be defined in Section 3, (3.8). Here we note that 
«3 = to 3 and 0:4 = to 4 — 2. The summation on the right-hand side of (2.9) is taken over 
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a finite set of non-negative integer pairs (p, m). The coefficients c PiTn can be calculated 
explicitly. For the cases k — 1, 2 we have 

BAz) = a 3 - 



and 



'1/z-z 

(l) = ( Q '*- Q '3i7r^ + a K(i7^ + (i7r-^)- {2M) 



Note that 



5i(G,J^)) = - 7 |^GjJ^) = -«3|_ 2 2 ^4^2(^/2K(a;)), z G C + . (2.11) 
If as — 0, then 

B 2 {G^{z)) = -j=^Gl B W zeC+. (2.12) 

Now we can formulate a counterpart of Edgeworth expansion in the Free CLT. We 
obtain this counterpart from the following results in which we establish an approximation 
of the measures fi n by the free Meixner measures. Consider the three-parameter family of 
probability measures {fi a ,b,d :oGl,6< l,d< 1} with the reciprocal Cauchy transform 

q 1 =a+ l -((l + b){z-a) + > /(i-6)2( g _ a )2_4( 1 _ d )^ zeC , (2.13) 

f^a,b,d V'/ 

which we will call the free centered Meixner measures (i.e. with mean zero). In this 
formula we choose the branch of the square root determined by the condition Qz > 
implies Q(l/G flabd (z)) > 0. These measures are counterparts of the classical measures 
discovered by Meixner [34]. The free Meixner type measures occurred in many places in 
the literature, see [3], [17], [18], [19], [29], [33], [37]. 
Assume that m 4 < oo, mi = 0, m 2 = 1 and denote 

m 3 , m 4 — ml — 1 , m 4 — m\ TNT 
a n :=^, b n :=^ 3 - , d n := -± 3 -, n G N. (2.14) 

In the sequel we will use the free Meixner measures of the form /io,o,o — w , Hanflfi if 
03 < oo, mi = 0, m 2 = 1 and [ia n ,b n ,d n if rri4 < oo, mi = 0, m 2 = 1 and n > m 4 . 

Recall that a probability measure n is ffl-infinitely divisible if for every n G N there 
exists u n G M. such that \i = u n ffl u n ffl • • • ffl z/ n (n times). 

Using the results of Saitoh and Yoshida [37], we will show in Section 4 that under 
the assumptions (3^ < oo and n > m| the free Meixner measure // a „,o,o is absolute contin- 
uous with a density of the form (4.1), where a = a n , b — 0, d — 0, and /i a „,o,o is ffl-infinitely 
divisible. Under the assumptions m 4 < oo and n > 3m 4 the free Meixner measure fi an ,b n ,d n 
is absolute continuous with a density of the form (4.1), where a = a n ,b = b n ,d = d n , and 
Va n ,b n ,dn is ffl-infinitely divisible. 
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We now introduce some further notations. Assume that f3 q < oo for some q > 2. 
Introduce the Lyapunov fractions 

L qn := , g 9V9 and let p q (fJ,,t) := / \u\ q /i(du), t > 0. (2.15) 
n (q 2)/2 J\u\>t 

Write 

gi := min{g, 3}, q 2 : = min{g, 4}, q 3 : = min{g, 5}. 
Then denote, for n E N, 

^ s (n):= inf ^(e), where ggM (e) := £ * +2 ~^ + Es^M^L ^ (2 .l 6 ) 

0<e<10- 1 /2 

provided that /3 q < oo, q > s + 1, for s = 1,2,3, respectively. It is easy to see that 
< Vqs( n ) < 10 1+s / 2 + lfors+l<g s <s + 2 and rj qs {n) — > monotonically as n — > oo 
if s + 1 < q s < s + 2, and i] qs (n) > 1, n e N, if q s — s + 2. 

By agreement the symbols c, Ci, c 2 , . . . and c(/i), Ci(/i), c 2 (/x), . . . shall denote absolute 
positive constants and positive constants depending on ji only, respectively. By c and 
c(n) we denote generic constants in different (or even in the same) formulae. The symbols 
ci,C2, ■ ■ ■ and ci(/i), C2(/i), . . . are applied for explicit constants. 

Theorem 2.1. Assume that Xj, j = 1,..., are free, f3 q < oo with some q > 2 and 
mi = 0, m 2 = 1. Then, for n£N, 

sup |F„ W - „(*)| < c + ^ < 2 < 3 (2.17) 

x-elR [^3n, y P g < oo, q > 3. 

In the case m 2 < oo, Theorem 2.1 yields a type of Free CLT with the error bound 
sup \F n (x) — w(x)\ < c(i] ql (n) + n' 1 j , fiGN. 



Since i] q i(n) < 10 3//2 + 1, n E N, in the case /3 (? <oo,2<g<3, we obtain from (2.17) 
the complete analogue of the Berry-Esseen inequality as well. 

Corollary 2.2. Assume that Xj, j = 1,..., are free, (3 q < oo with 2 < q < 3 and 
mi = 0, m 2 = 1. Then, for n E N, 

sup \F n (x) - w(x)\ < cL qn . (2.18) 

x-eR 

In the case f3 3 < oo the inequality (2.18) has the form 

sup \F n (x) - w(x)\ < cL 3n , n E N. (2.19) 
x-eR 

The upper bound (2.19) sharpens previous results obtained by the authors [22] and V. 
Kargin [26]. 

Theorem 2.1 and Corollary 2.2 are free analogues of Esseen's inequality in classical 
probability theory (see [36], p. 112-120). 
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Theorem 2.3. Assume that Xj, j = 1, . . . , are free, f3 q < oo with some q > 3 and 
mi = 0, m 2 = 1. TTien, /or n£N, 

sup \F n(x) - ^,„,„((-oc,,))| < c + LL * ft<co,3S,<4 

Corollary 2.4. Under the assumptions of Theorem 2.3 £/ie following expansion holds 

F n {x) = w(x) - -a n U 2 (x/2)p w (x) + p nl (x), xeR, (2.21) 

where the remainder term p n i(x) admits the bound, for x G R, n G N, 

U f-rU <r J%(«)V + 4 + K| 3/2 «/ & < oo, 3 < q < 4 
[L 4n + lanl^ 2 i/ /3, < oo, q > 4. 

Note that in the case f}$ < oo the estimate (2.22) yields the bound 

ki(z)| < c(n q2 (n) + L 3n + |a n | 1/2 )L 3n , (2.23) 

where r] q2 (n) — > as n — > oo, and we obtain an analogue of Edgeworth expansion. 

Since rj q2 (n) < 101, 3 < q < 4,n G N, the results (2.21) and (2.22) again yield the free 
Berry-Esseen inequality (2.19) as well. 

In addition we obtain from Theorem 2.3 the following bounds. 

Corollary 2.5. Under the assumptions of Theorem 2.3 

sup |F n (a;) - /i anA o((-oo,a;))| < cL qn for n G N if P q < oo, 3 < q < 4. (2.24) 

Before formulating the next result, denote by q n , n > m 4 , a signed measure with the 
density 

P^( x ) '■= ( e l( x ~ a nf ~ l)p w {e n {x - On)), xER, (2.25) 
where e n := (l — b n )/y/l — d n . Denote by K n , n > m 4 , the signed measure K n := fi an ,b n ,d n + 
^q n . It is easy to see from results of Section 4 that n n is a probability measure for n > m 4 / c 
with some sufficiently small c. 

Theorem 2.6. Assume that Xj, j — 1, . . . , are free random variables, that (3 q < oo with 
some q > 4 and that m\ — 0, m 2 = 1. Then, for n > m 4; 

i jji / \ // ui . )Vq3(n)L qn + L 4 ^ 2 i/ /5, < oo, 4 < q < 5 
sup \F n (x) - K n ((-oo,x))\ < c< q y (2.26) 

xeR [L 5n if fJ q < oo, q > 5. 

Corollary 2.7. Assume that the assumptions of Theorem 2.6 are satisfied. Then 
F n (x + a n ) = w(x) 

+ ( - ^(x/2) + ^(3 - U 2 (x/2)) - bn ~ a \~ 1/n U 3 (x/2))p w (x) + p n2 (x), 

(2.27) 
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for all real x, where 

M*)\<ch 3(n)Lv ' + Ll " ^' <00 ' 4 ^ <5 for I£ l neN. (228) 
[L 5n if p q < oo, q > 5 

If m 3 = this formula has the following simple form 

F n (x) = w(x) - ^^-U 3 {x/2)p w (x) + p n3 (x), (2.29) 
where p n3 (x) admits the bound (2.28). 

If m 3 7^ 0, we obtain from (2.27) the following expansion for F n (x): 

F n (x) = w(x) - l -a n U 2 {x/2)p w {x) + - K ~ ^ ~ 1/n U 3 (x/2))p w (x) 

+ Qi{x, a n ) + Q 2 (x, a n , b n , l/n) + p ni (x), x e R, 

where 

Qi(x, a n ) = w(x - On) - w(x) + a n p w (x) + y (3 - U 2 (x/2))(p w (x - a n ) - p w (x)), 

Q 2 (x,a n ,b n ,l/n) = (^tfi(s/2) - ^ ~ °\~ 1/n U 3 (x/2)) (p w (x - a n ) - p w {x)) 

and the function pn^x) admits the bound (2.28). The function Qi(x,a n ) is a function of 
bounded variation and it is not difficult to verify that 

-M 3/2 < sup|Qi(a;,a n )| < c\a n f 2 and -\a n \ 3/2 < \\Qi{x,a n )\\ TV < c\a n \ 3/2 , (2.30) 

with some c > 1. This means that Qi(x,a n ) is actually of order n~ 3 / 4 . We shall see 
that Qi(x, a n ) can not be cast by Taylor expansion around x into an expansion in powers 
of n _1//2 like (2.27) in terms of p w (x) and the Chebyshev polynomials which is continu- 
ous up to the boundary ±2 with finite total variation. Moreover \Q 2 (x,a n ,b n ,l/n)\ < 

Indeed from the formal expansion of p n in (2.8) and (2.11) it follows that the first 
two summands on the right-hand side of (2.8) are Cauchy transforms of the finite signed 
measure on the right-hand side of (2.21). Moreover, in the case m 3 = the first three 
summands on the right-hand side of (2.8) are Cauchy transforms of the finite signed 
measure on the right-hand side of (2.29). But in the case m 3 ^ the third summand on 
the right-hand side of (2.8) can not be a Cauchy transform of a signed measure ( which 
is finite on every bounded interval and f R \((du)\/(l + \u\) < oo. This will be proved in 
Section 4. Therefore, taking into account the formal expansion (2.8), we can not expect 
an expansion of type (2.27) for the function F n {x) without shift. 

Remark 2.8. The methods used in the proof of Theorems 2.3 and 2.6 still do not yield a 
free analogue of Edgeworth asymptotic expansions under the assumption f3 q < oo, q > 5, 
with a remainder term of order 0(n _3//2_7 ) with 7 > 0. This problem remains open. 
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Remark 2.9. It is known, see for example [10], [21], that there is semigroup ji t G -M, t > 1, 
such that <t>^ t {z) = t^^z), where 4>^ t {z) are Voiculesku transforms of the probability 
measures fi t - For the definition of Voiculesku's transform, see in Section 3. As before let 
mi = and m 2 = 1. Define a probability measure fi t in the following way: fi t ((—oo,x)) = 
fj,((—oo,xy/t)), iel Theorems 2.1, 2.3, 2.6 and their Corollaries remain valid for fi t if 
the integers n are replaced by t > 1. One can prove these results exactly by the same 
proof. 

It was proved in [5] that if the distribution \x of X x is not a Dirac measure, then 
F n (x) is Lebesgue absolutely continuous when n > ci(/i) is sufficiently large. Denote by 
p n (x) the density of F n (x). If /i has a compact support, Voiculescu's result [40] shows 
that the support of F n (x) is contained in the interval [—2 — 4=, 2 + -^=] for n > 1, where 
L := sup{|x| : x G supp(fj,)} (see [27] as well). Our method allows to obtain an asymptotic 
expansion for p n (x) in this case. 

Theorem 2.10. Assume that ji has compact support and mi = 0, m 2 = 1. Then, for 
n > c i(/-0> Pn(x) is a continuous function such that p n (x) < 2, x G R, and 

p n (x+a n )= (l + \d n -al---a n x-(b n -al--)x 2 )p w (e n x)-\ — — — *^E== (2.31) 
V 2 n V nJ J n 3 ' 2 ^4 - (e n x) 2 

for x G [—j- + h, j — h], where h = ^yr- Moreover, as a simple consequence of (2.31), 
the following inequality holds 



J 



Pn (x + a n ) dx < (2.32) 



Here and in the sequel we denote by 9 a real- valued quantity such that \9\ < 1. 
We see that the remainder term in (2.31) has the order — k= on the interval [— — + 
Mn A _ Mjl] w ith any M n such that M n -> oo and ^ -)• as n -)• oo. 
Now we shall state some consequences of Theorem 2.10. 

Corollary 2.11. Let fi n be the distribution of Y n from (2.1) with bounded identically 
distributed free summands X 1 , . . . , X n such that mi = and m 2 = 1. //m 3 7^ 0, then 

p n {x)-p w {x)\dx= i ^ + 9(c\a n f 2 + '^), n>cM. (2.33) 



/ 

Jr 



7T\/n n 

If m 3 = 0, then 

2|m 4 -2| c(/i) 

h (7— 

7rn n 



/ 



Pn(g)-P«,(g)l<fa= _„ +^^^ n>ci0i). (2.34) 



Recall that, if the random variable X has density /, then the classical entropy of a 
distribution of X is defined as h(X) = — J R f(x) log f(x) dx, provided the positive part 
of the integral is finite. Thus we have h(X) G [—00, 00). 



10 



G. P. Chistyakov and F. Gotzc 



A much stronger statement than the classical CLT - the entropic central limit theorem 
- indicates that, if for some n , or equivalently, for all n > n , Y n from (2.1) have 
absolutely continuous distributions with finite entropies h(Y n ), then there is convergence 
of the entropies, h(Y n ) — > h(Y), as n — > oo, where Y is a standard Gaussian random 
variable. This theorem is due to Barron [4]. Recently Bobkov, Chistyakov and Gotze [16] 
found the rate of convergence in the classical entropic CLT. 

Let v be a probability measure on R. The quantity 

x( v ) = / / log 1^ -y\v{dx)v(dy) + ^ + ]- log27r, 

called free entropy of u, was introduced by Voiculescu in [43]. Free entropy x behaves 
like the classical entropy h. In particular, the free entropy is maximized by the standard 
semicircle measure fj, w with the value xi^w) = |log27re among all probability measures 
with variance one, see [25], [44]. Wang [46] has proved the free analogue of Barron's result. 
We give the rate of convergence in the free CLT for bounded free random variables. 

Corollary 2.12. Let fi n be the distribution of Y n from (2.1) with bounded identically 
distributed free summands Xi, . . . , X n such that mi = and m<i = 1. For n > ci(/i), 

r r 3 1 m 2 c(a) 

X(l*n) =ii log [a; - y\p n (x)pn(y) dxdy + - + - log 2tt = x(fJ> w ) ~ + d ^H- 
J Jrxr 4 z on n > 

Suppose that a measure v has a density p in L 3 (R). Then, following Voiculescu [43], 
the free Fisher information of v is 

,2 



$(z/) = ^— [ p(xf dx. 
3 Jr 



It is well-known that $(w) = 1. Moreover, the free Cramer-Rao inequality shown in 

[43] says that J R (x — J R uv{du) s j v[dx) > 1, and equality holds if and only if v 

is a measure with a semicircle distribution function. We obtain the following result for 
bounded free random variables. 

Corollary 2.13. Let fi n be the distribution of Y n from (2.1) with bounded identically 
distributed free summands Xi, . . . , X n such that mi = and m 2 = 1. For n > ci(/i), 

*W = ^ I Pn{xfdx = ^ w ) + r ^ + C ^l. (2.35) 



3. Auxiliary results 



We need results about some classes of analytic functions (see [1], Section 3, and [2], 
Section 6, §59). 
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The class Af (Nevanlinna, R.) is the class of analytic functions f(z) : C + — > {z : ^sz > 
0}. For such functions there is the integral representation 

f 1 ~\~ uz f / I u \ 

f( z \ = a + oz + / r(du) = a + bz + / (1 + w 2 ) r(<iw), z G C + , 

J R u- z J R \u-z 1 + U 2 J 

(3.1) 

where i) > 0, a G I, and r is a nonnegative finite measure. Moreover, a = $lf(i) and 
r(R) = — 6. From this formula it follows that 

f(z) = (b + o(l))z (3.2) 

for z G C + such that stays bounded as \z\ tends to infinity (in other words 

z — > oo nontangentially to R). Hence if 6 7^ 0, then / has a right inverse defined on 
the region 

r Q)/3 := {2 G C + : < aSz, 3z > /?} 
for any a > and some positive /3 = f3(f,a). 
A function / G Af admits the representation 

, 6 C+, (3.3) 

where a is a finite nonnegative measure, if and only if sup y>1 \yf(iy)\ < 00. 

For ji G A4, consider its Cauchy transform G^(z) (see (2.7)). The measure /i can be 
recovered from G^(z) as the weak limit of the measures 

fi y (dx) = ^sG^(x + iy) dx, i£l, y > 0, 

as y j. 0. If the function QG^z) is continuous at rr G R, then the probability distribution 
function D^(t) = fJ>((— 00, £)) is differentiable at a; and its derivative is given by 

D'^x) = -SG^/ir. (3.4) 

This inversion formula allows to extract the density function of the measure \i from its 
Cauchy transform. 

Following Maassen [30] and Bercovici and Voiculescu [9], we shall consider in the fol- 
lowing the reciprocal Cauchy transform 

F,W = a^y (3.5) 

The corresponding class of reciprocal Cauchy transforms of all \i G A4 will be denoted by 
J 7 . This class coincides with the subclass of Nevanlinna functions / for which f(z)/z — > 1 
as z — > 00 nontangentially to R. Indeed, reciprocal Cauchy transforms of probability 
measures have obviously such property. Let / G Af and f(z)/z — > 1 as z — > 00 nontan- 
gentially to R. Then, by (3.2), / admits the representation (3.1) with 6=1. By (3.2) 
and (3.3), —l/f(z) admits the representation (3.3) with a G A4. 
The functions / of the class J 7 satisfy the inequality 

%f(z) >%z, z G C + . (3.6) 
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The function <j>n(z) = *\z) — z is called the Voiculescu transform of fi and <p^{z) 
is an analytic function on Y a ^ with the property 9^(2) < for z E F a ,/3, where <p^(z) 
is defined. On the domain T a ^, where the functions <p^ 2 (z), and Ml ffl M2 (<2) are 

defined, we have 

<P^m^(z) = <j>n(z) + <f>ua(z). (3.7) 

This relation for the distribution /j,i ffl ji 2 of X + Y, where X and Y are free random 
variables, is due to Voiculescu [40] for the case of compactly supported measures. The re- 
sult was extended by Maassen [30] to measures with finite variance; the general case was 
proved by Bercovici and Voiculescu [9]. 

Assume that (3 k < oo for some k EN. Then 

n I \ 1 m l m k ( 1 \ „ 

G mW = ; + ^ + --- + ^h+«(^^TtJ ) z-Kx>, zET aA . 
It follows from this relation (see for example [27]) that 

a 2 a k ( 1 \ „ 

^(z) = «! + — + ••• + + °( -jfc=T J' z ~* 00 ' 2:6 r °' 1 - ^ 3 - 8 ) 

^ Z \ Z / 

We call the coefficients a m , m = 1, . . . , k, the free cumulants of the probability measure 
fi. It is easy to see that a± — mi, a 2 — rn 2 — rn\, a^{ji) = m 3 — 3mim 2 + 2m\ . In the case 
mi = and m 2 = 1 we have ot\ = 0, 0:2 = 1, «3 = m 3 and 0:4 = m 4 — 2. 

If /i G M. has moments of any order, that is (3 k < 00 for any A; G N, then there exist 
cumulants a m , m — 1, . . . , and we can consider the formal power series 



m=l 

In addition M (z) satisfies (3.8) for any fixed k E N. If fi has a bounded support, (f>n(z) 
is an analytic function on the domain \z\ > R with some R > and the series (3.9) 
converges absolutely and uniformly for such z. 

Voiculescu [43] showed for compactly supported probability measures that there exist 
unique functions Zi,Z 2 E T such that G^m^z) = G^Z^z)) = G^ 2 (Z 2 (z)) for all 
z E C + . Using Speicher's combinatorial approach [38] to freeness, Biane [15] proved this 
result in the general case. 

Chistyakov and Gotze [21], Bercovici and Belinschi [6], Belinschi [7], proved, using 
complex analytic methods, that there exist unique functions Z\(z) and Z 2 (z) in the class 
T such that, for z E C + , 

z = Zi(z) + Z 2 (z) - F^Ziiz)) and F^Z^z)) = F, 2 {Z 2 {z)). (3.10) 

The function F lll (Zi(z)) belongs again to the class J 7 and there exists \x E A4 such that 
F Ml (Zi(z)) = Fn(z), where F^z) = 1/G^z) and G^(z) is the Cauchy transform as in 
(2.7). We can define the additive free convolution in the following way /ii ffl /i 2 := fi. 
The measure fi depends on fi\ and fi 2 only. The relation (3.7) follows immediately from 
(3.10) and we see that this definition coincides with the Voiculescu, Bercovici, Maassen 
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definition. Hence we have the equivalence of a "characteristic function" approach and 
a probabilistic approach to the definition of the additive free convolution. 

Specializing to fj,± — /i 2 — • • • — — A* write fj,i ffl • • • ffl fj, n — fi nm . The relation (3.10) 
admits the following consequence (see for example [21]). 

Proposition 3.1. Let /i G M.. There exists a unique function Z G T such that 

z = nZ(z)-(n-l)F fM (Z(z)), z G C + , (3.11) 
and F^ n m(z) = F /Ji (Z(z)). 

The next lemma was proved in [21]. 
Lemma 3.2. Let g : C + — > C~ be analytic with 

liminfM!^ = . (3.12) 

y^+oo y 

Then the function f : C + — > C defined via z h- >■ z + g(z) takes every value in C + precisely 
once. The inverse : C + — > C + thus defined is in the class J 7 . 

This lemma generalizes a result of Maassen [30] (see Lemma 2.3). Maassen proved 
Lemma 3.2 under the additional restriction \g(z)\ < c{g)/^sz for z G C + , where c(g) is 
a constant depending on g. 

Using the representation (3.1) for F^(z) we obtain 

F„(z) = z + XF»{i)+ t (l + Uz)T(iU \ :E C», (3.13) 

Jr u - z 

where r is a nonnegative measure such that r(R) = StF^i) — 1. Denote z = x + iy, where 
x, y G R. We see that, for ^sz > 0, 

3 (nz - (n - l)F^z)) = y(l - (n - l)I^(x,y)\ where I^(x, y) := f ^ + u . 
V /V / J R (u - x) 2 + y 2 

For every real fixed x, consider the equation 

y(l -{n- l)I^x, y)) = 0, y > 0. (3.14) 

Since y >->■ I^x, y), y > 0, is positive and monotone, and decreases to as y — > oo, it is 
clear that the equation (3.14) has at most one positive solution. If such a solution exists, 
denote it by y n (x). Note that (3.14) does not have a solution y > for any given x G R 
if and only if I^x, 0) < \/{n - 1). Consider the set S:={ieR: I^x, 0) < l/(n - 1)}. 
We put y n (x) = for x G 5. By Fatou's lemma, / M (x ,0) < liminf^^ I^x, 0) for 
any given rr G R, hence the set 5 is closed. Therefore R \ S is the union of finitely 
or countably many intervals (xk,Xk+i), %k < x k+i- The function y n (x) is continuous on 
the interval (xk,Xk+i)- Since the set {z G C + : nQz — (n — 1)9^(2:) > 0} is open, we 
see that y n (x) — > if x j. Xk and x t ^fc+i- Hence the curve 7„ given by the equation 
z = x + iy n (x), x G R, is continuous and simple. 

Consider the open domain D n := = x + iy, x, y G R : y > y n (x)}. 
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Lemma 3.3. Let Z E J 7 be the solution of the equation (3.11). The function Z(z) maps 
C + conformally onto D n . Moreover the function Z(z), z E C + ; is continuous up to the real 
axis and it establishes a homeomorphism between the real axis and the curve j n . 

Proof. We obtain from (3.11) the formula 

Z { - l \z)=nz-{n-l)F^z) (3.15) 

for z E r Q)/ g with some a, (3 > 0. By this formula we may continue the function Z^~ x \z) 
as an analytic function to C + . Using the representation (3.13) for the function F^z), 
we note that Z^~^(z) = z + g(z), z E C + , where g(z) is analytic on C + and satisfies 
the assumptions of Lemma 3.2. By Lemma 3.2, we conclude that the function Z^~ x \z) 
takes every value in C + precisely once. Moreover, as it is easy to see, Z( _1 ^(D n ) = C + . 
The inverse Z(z) gives us the conformal mapping of C + onto D n . By the well-known 
results of the theory of analytic functions (see [32]), Z(z) is continuous up to the real axis 
and it establishes a homeomorphism between the real axis and the curve j n . □ 

Lemma 3.4. Let n be a probability measure such that mi = 0,m2 = 1. Assume that 
p 2 (A*, \f{n — l)/8) < 1/10 for some positive integer n > 10 3 . Then the following inequality 
holds 

\Z(z)\ > y/(n- l)/8, zEC + , (3.16) 
where Z E T is the solution of the equation (3.11). 

Proof. Write 

^ , n 1 r ( z ) ^,4. / \ f u 2 u(du) 

GJz) = - + -V, z E C where r(z) : = / PV . 3.17 
z z 2 J R z-u 

It is obvious that \r(z)\ < l/^sz, z E C + . Rewrite (3.14) in the form 

9 ,(l + ( n -l)(l-^^))=0. 

Let us show that 

y n (x) > y/(n-l)/8 for \x\ < y/(n - l)/8. (3.18) 
In order to prove this inequality we shall establish that 

( "- 1) ( 1 -^ 9 G^)) < - 1 (3 ' 19) 
for \z\ = \yj{n — 1) and \^tz\ < y/(n — l)/8. Indeed, since the function 

-J M (^,^) = (n-l)(l-^S-^), Zz>0, 

\ ^SZ Lr^{Z)/ 

is negative and monotone, and increases to as ^sz — > oo, (3.18) follows from (3.19). 
We have, for the z considered above, 

1 (I r( z y 



Gn(z) \z z 



(- + ArJ = *-r(*) + n(*), 
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where r-y(z) admits the upper bound |ri(z)| < 2(\z\(^sz) 2 y 1 < 32/ (n — l) 3 / 2 . 
Using the previous formula, we easily obtain the relation, for the same z, 

-/,(*, 9,) = (n - + „(,) = -(„ - 1) I ( „_^ff (sb) , + r,W, (3.20) 

where r 2 (z) admits the upper bound |r 2 (z)| < 32VS/(n - 1) < 1/6. Hence we have, for 
the same z, 

< -c-i)/^ +r ; M 

< (1 - P 2 (p, + r 2 (*) < " P 2 (/i, V(n - l)/8)) + r 2 (z) 

6 1 

< -5 + 6 < " 1 
and (3.19) is proved. 

The assertion of the lemma follows immediately from (3.18). □ 

Denote by A(/c', k") the Kolmogorov distance between the finite signed measures k' 
and k" such that k'((—oo,x)) — >■ and «"((— oo, #)) — >■ as x — >■ — oo, i.e., 

A(k',k") := sup |k/((-oo,x)) - «"((-oo,a:))|. 

We need the following result of Bercovici- Voiculescu [9] . 
Proposition 3.5. If fj,,fj/,i/ and v' are probability measures, then 

A(n mv,n'm v) < A(/i, //) + A(v, is'). 
In addition the following proposition holds (see [36], p. 139). 
Proposition 3.6. If 3 < m < k, then the Lyapunov fractions L mn and Lk n satisfy the 

-l/(fc-2) 

4. Properties of free Meixner measures 



inequality: Llln 1 ^ < L 1 ^ 1 ' 



Saitoh and Yoshida [37] have proved that the absolutely continuous part of the free 
Meixner measure fj, a ,b,d, a e R, 6 < 1, d < 1, is given by 



y/4(l-rf)-(l-fe) 2 (x-a) 2 

2tt/(x) ' l4 ' ij 



when a - 2y/l - d/(l - b) < x < a + 2y/l - 3/(1 - b), where 

f(x) := bx 2 + a(l - b)x + 1 - d; 

the measure may have a discrete part /ip in the following cases: 
1. if /(x) has two real roots yi 7^ y 2 , then 

Hd := Ai<J yi + A 2 <5 W , (4.2) 
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where 

1 



A, 



^a\l-bf-Ab{l-d)\ \yj\ 
2. if 6 = and a ^ 0, then 

l-d\ „ , l-d 



:= (l — ) 5 y , where y := . (4.4) 

Recall that 5 y with y e K. is a Dirac measure concentrated at the point y. 

Saitoh and Yoshida proved as well that for < b < 1 the (centered) free Meixner 
measure fi a ,b,d is ffl-infinitely divisible. Note (see Bozejko and Bryc [17]) that n a ,b,d — l^w 
if a = b = d = 0; (J, a ,b,d is the free Poisson type measure, which is also known as Marchenko- 
Pastur measure [31], if b — d — and a ^ 0, and fi a ,o,d with a ^ 0, d ^ is the shifted 
free Poisson type measure ; yt, a ,b,d is the free Pascal (negative binomial) type measure 
if b > and a 2 (l — b) 2 > 46(1 — d); fi a ,b,d is the free gamma type measure if b > 
and a 2 (l — b) 2 = 46(1 — d); n a ,b,d is the pure free Meixner type measure if 6 > and 
a 2 (l -6) 2 < 46(1 -d). 

Now assume that m 4 < oo, mi = 0, m 2 = 1 and n > 3m 4 . By the well-known moment 
inequality 



1 mi m 2 
mi m 2 W3 
m 2 m3 



> 



(see [1]), we conclude that m 4 — 1 — m| > 0. Therefore the lower bounds b n > and 
d n > hold. In addition we have \a n \ < 1/V3 : b n < 1/3 and d n < 1/3. Consider 
the measures fi an ,b n ,d n - These measures may be the free Pascal, the free gamma and the 
pure free Meixner type measures. 

Let 6 n > 0. Note that the polynomial f n (x) = b n x 2 + a n (l — b n )x + 1 — d n has two real 
roots yi n and y 2n in the case a 2 (1 — 6 n ) 2 — 46 n (l — d n ) > and these roots have the same 
sign. By the relation 

1 1 _ H(1-6 TO ) < yg < 1 
1 2/i I | J/2 1 1-4 2 
one can deduce the inequalities |?/j„| > 1, j = 1,2. Using (4.2) and (4.3) we see that 
the discrete part of fi a „,b n ,d n is equal to zero. Let b n = and a n ^ 0. We see from (4.4) 
that in this case the discrete part of fj. an ,b n ,d n is equal to zero as well. 

Thus, in the considered case it follows from Saitoh and Yoshida's results that the prob- 
ability measure Ha n ,b n ,d n is ffl-infinitely divisible and it is absolutely continuous with a den- 
sity of the form (4.1) where a = a n ,b = b n ,d = d n . 

Assume that (3% < oo,mi = 0,m 2 = 1 and n > m 2 , i.e., \a n \ < 1. In this case 
the probability measure /i a „,o,o is absolute continuous with a density of the form (4.1) 
where a = a n , 6 = 0, d = 0. In addition, by Saitoh and Yoshida's results, (J, a „,o,o is 
ffl-infinitely divisible. 
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5. Formal Asymptotic expansion in the Free CLT 
In this section we deduce formula (2.8). 

By Proposition 3.1, there exists Z(z) G T such that (3.11) holds, and F fin m(z) = 
F^Z(z)). Hence F^(z) = F^S n (z))/^i, z G C+ where S n (z) := Z{^iz)/^. 
Using the Voiculescu transform (see (3.7), this relation implies that 

S n (z) = F^z) + MV^F^z))/^ 
for z G r Q)i a with some a, (3 > 0. On the other hand we conclude from (3.11) that 

7 71 — 1 

S n (z) = - + F, n (z), zeC + . 

n n 

The last two equations give us 

F^z) + yM^F^z)) = z, zeY afi . (5.1) 

Consider the function f(z) := z + y/nfi^y/nz), z G r Q ( g/ with some (3' > (3, and define 
the function 

:=!(/(*) + V/ 2 (*)- 4), ^r a/) (5.2) 

where we choose the branch of the square root by the condition ^sg(z) > for z G T a ^>. 
It is easy to see that g(z) = z(l + o(l)) as z — > 00 nontangentially to R. In addition, by 
(5.1), g(z) satisfies the relation 

We deduce from (5.3) that 

9(FM) = l(z + V^Z)=F ltm {z) 1 zeTajy. 

Recall that we denote by fi w the semicircle measure. 

Since the function g(z) has a right inverse g^~ 1 \z) in Y a ^» with some (3" > (3', we have 

F lin (z)=g(- 1 \F lia (z)), ^V, where /T > (3". (5.4) 

Let /i G M. such that all moments of \i exist. In addition let m\ — and m 2 = 1. 
Consider the formal power series in z 

00 

vW„(vM := ^ (fc -iv2 fc , (5-5) 
fc=i 

where /c = 1,2,..., are free cumulants of the measure \i and the formal power series 
of g: 

00 

*(*) = *+£5- (5 - 6) 

k=0 
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In our case a± = 0,a 2 = l,c*3 = tuz and 0:4 = m 4 — 2. Using (5.3) and (5.5), (5.6) we 
obtain the following relation for the considered formal power series 



k=0 k=0 k=0 k=l 

It follows from this relation that ao = 0, a\ = and 

fe-3 

\fc-l„fc-l a k+i 



a k - ak-2 + ^a s a fc - s -3 h (~l) fc ^0 1 = ( fc _"t)/2 ' fc = 2,3,.... (5.7) 

s=0 ^ 

We have from (5.7) the relations a 2 = otz/y/n, a 3 = a 4 /n. In addition we obtain from 
(5.7) by induction that 

°- = ^ + 0OA and Bw= «t- ('- 1 )(«- 2 ) M + (-^) (5.8) 

as n 00 for s = 2, ... . 

Now consider the formal power series for the right inverse g ly ~ 1 \z) 

fc=0 



2+> h> r. — Z. 

m=0 



Rewrite the relation g(g( ^(z)) = z in the form 

00 

Y.—, — 



Using the formula 



we finally get 



OO 7 OO OO /7-i,\°°7 



m=0 fc=2 s=0 v 7 m=0 

We obtain from this equality that b = b\ = 0, b 2 = —a 2 and 

m ~ 1 m ~ k /ft _ 1 + s \ 
& m + a m + ^a fe ^(-lW J 6mi • • • 6m. = 0, m = 3, . . . . 

k=2 s=l ^ ' miH \-m s =m— k— s 

(5.9) 

Moreover it is easy to deduce from (5.8) and (5.9) that 

^=-«-+°y 7i )=-^! +0 (j 7J ) (5.10) 
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and 



m— 2 



t>2m-l — —a2m~l + 2 Sa2s^2m-2-2s + — 3/2) 

s=l 

= _g._ (m-2)(m + l) 0j + o /_M m= (5U) 
ra 2 n \n d ' 2 / 

for m = 2, . . . . It remains to note that 

and we can write the formal power series in l/y/n 

1 J B k (l/ Z ) 

9 { - 1] ( z ) z t^i nk/2 ' 
Taking into account the relations (5.10) and (5.11), we easily conclude that 

°° 1 11 
B 1 (l/z) = a 3 ~^ = "3^3 • j^j 



m=2 

and 



„ . . . ^ 1 9 / (m — 2)(m + 1) 1 1 / ^ 1 V 

m=2 rrt=2 m=l 



Since 

E 7/t ± / \ 1 

r 2m+l 9 V Z / 7 2m 



v m 1,^ 1 V 1/ 1 V 1 1 



z 2m+i 2 V^z 2m / 2V^ 2 (z 2 -l)/ ^ 3 (z 2 -l) ^(^ 2 -l) 2 ' 

m=2 m=2 v ' \ / \ / 



E77T 1 / v - m\' 1 1 

^Tr^ _ 2V^^J = z 3 (z 2 - 1) + ^(^ 2 - l) 2 + (z 2 - l) 3 ' 



2z 



m=2 m=2 

we finally obtain 



In view of (5.7) and (5.9), we see as well that B k (z) are functions of the form 

1 1 



B k (l/z) = 



'zp(z- l/z) m 

with real coefficients c Pi?n which depend on the free cumulants 0:3, ... , a^+2 and do not 
depend on n. The summation is carried out over a finite set of non-negative integer pairs 
(j), m). The coefficients c PtTn can be calculated explicitly in the way described above for 
the coefficients c P:fn of the functions B\(\/z) and B 2 (l/z). 
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Hence we deduce from (5.4) the formal expansion 

GM = G^z) + £ Bk{ % {z)) - (5-12) 



n k/2 

k=l 



Using integration by parts, it is not difficult to verify that 

[ —d(\u 2 (x/2) Pw {x)), zeC + . (5.13) 
J_ 2 z-x V3 / 



= -«3 

On the other hand we see that if 0:3 7^ then the function B 2 {G^ w {z)) is not the Cauchy 
transform of some signed measure. Indeed, it is easy to see, using direct calculations, that 

a 2 

B 2 (G» w (z)) = {z2 _\ )3/2 +g(z), ze£\ (5.14) 

where the function g(z) is analytic on C + and there exists finite limit, for every —00 < 
ti < t 2 < +00, 

lim / ^sg(x + iy) dx. (5.15) 



ti 



In addition we note that 

lim 



j 3/2 



Assume now that B 2 {G^ w {z)) is a Cauchy transform of a real- valued function u(x) of 
bounded variation on every bounded interval and such that 

00 

\d(jj(x)\ 

< 00. 



If 



+ Fl 

Then, by Stieltjes-Perron's inversion formula [1], we have 
ou(t 2 + 0) - ou(t 2 - 0) ufa + 0) - w(ti - 0) 1 



= lim- [ 2 QB 2 (G flw (x + ty))dx. (5.17) 



Assuming in (5.17) ti := 3/2 and i 2 : = 2, and taking into account (5.14)-(5.16), we arrive 
at a contradiction. 
If a 3 (/i) = 0, then 

B 2 (G, (z)) = * Gj (z) = ^f x4 7 4x2 + 2 

= -a 4 / d(^U 3 (x/2)p w (x)), zeC + . (5.18) 

J-2 Z — X \4: / 
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6. Edgeworth Expansion in the Free CLT (the case (3 q < oo, q > 3) 
In this section we prove Theorem 2.3 and Corollary 2.4. 

Proof of Theorem 2.3. Recall that we denote by /i n the distribution of Y n in (2.1) 
for the free random variables Xj. Our first step is to reduce the problem to the case of 
bounded free random variables. 

6.1. Passage to measures with bounded supports. Let n G Af. Let e n G (0, 10 -1 / 2 ] 
be a point at which infimum of the function g qn 2{£) from (2.16) is attained. This means 
that 

Without loss of generality we assume that 

i] q2 (n)L q2n + L 3n < d, (6.1) 

where C\ > is a sufficiently small absolute constant. By Lyapunov's inequality f3 3 > 
m^ 2 = 1, we obtain from (6.1) that n in this case has to be sufficiently large, i.e., 
n > q 2 /3f > c{ 2 . 

Consider free random variables X,Xi,X 2 , . . . with distribution fi = C(X) such that 
j2([—e n y/n,e n y/n\) = 1 and fi(B) = fi(B) for all Borel sets B C [—e n ^/n,e n ^/n\ \ {0}. 
Denote by fi n distribution of the random variable Y n := (X 1 + • — h X n )/y/n. In addition 
introduce random variables 

v* X ~ An v* ■ Xl ~ An v* ■ X<2 ~ Afl v* ■ 1 X * n 

A .= — ,A 1 .- — ,A 2 .- — and Y n — j= , 

<_/„ <_/„ <_/ n y n 

where 

2\ 1/2 



4i 



[ u fi(du) and C n := (l— j u 2 ji{du)—( f u ji{du)\ \ 

J\u\>£ n% /n V J\u\>£ n Jn V J \u\>£ n Jn 



'\u\>£ ny /n v J\u\>£ n% /n y J\u\>£ ny /n 

Denote by /i* and /i* the distributions of the random variables X* and Y*, respectively. 
We denote by m* k and rh k , k = 0, 1, ... , the moments and by (3%. and k — 0, 1, . . . , 
the absolute moments of the distributions fi* and ft,, respectively. It is obvious that m\ = 
and m\ = 1. Using (6.1) we note that 

14,1 < e-^-Vn-te-Wp^eny/n) < -^ Vq2 (n)L q2n (6.2) 

and 

< ^- - 1 < 2(p 2 (/i, e n ^n) + A 2 n ) < 3r] q2 {n)L q2n . (6.3) 
By (6.1)-(6.3), we obtain 

C-\e n ^+\A n \)< l -^i. (6.4) 
It follows from (6.4) that the support of /i* is contained in [— \y/n, \ \/n\- 
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By (6.1)-(6.3), we easily deduce as well that 

\m* 3 - m 3 1 < C; 3 \m 3 - m 3 \ + (C- 3 - l)\m 3 \ + C; 3 (?,\A n \m 2 + 3^|mi| + \A n \ 3 ) 

4 

< C~ 3 \m 3 - m 3 \ + 4\m 3 \r) q2 (n)L q2n + —rj q2 (n)L q2n 

V n 

< C- 3 e-^n-^- 3 ^ 2 Pq2 (fi,e n ,/E) + 4(jm 3 | + -£)r) q2 (n)L q2n < 2^r lq2 {n)L q2W 

(6.5) 

and, using similar arguments, 

4 

fit < C~ 3 f3 3 + —rj q2 (n)L q2n , m* < C~ 4 m 4 + 5L 3n rj q2 (n)L q2n . (6.6) 



n 



Let T be a random variable with distribution /i an ,o,o- Denote by /i an ,o,o the distribution 
of C n T + ^/nA n . 

By the triangle inequality, we have 

A(/j„, /i an ,o,o) < A(/i n , /i n ) + A(fl n , At an ,o,o) + A(/2 a „ A0 , /V,,o,o)- (6.7) 
First we establish with the help of Proposition 3.5 

A(// n ,// n ) < nA(/i,/i) < n/j({\u\ > e n y/n}) < e^ q2 n~^ 2)/2 p Q2 (ii, e n y/n) < r] q2 (n)L q2n . 

(6.8) 

Recalling the definition of /i a „,o,o (see (2.13) and (4.1), (4.4)), we note that p an ,o,o is 
an absolutely continuous measure with the support on [a n — 2,a n + 2] and its density has 
the form 



4 — (x — a n ) 2 /{2n{l + a n x)), when x G [a n — 2, a n + 2]. (6.9) 

This density does not exceed 1 on the set [a n — 2, a n + 2] and is equal to outside of this 
set, therefore we easily deduce the following upper bound, using (6.2) and (6.3), 

/ 1 \/nA \ 

A(/2 a „,o,o, Va n ,o,o) < ci^— - 1 + —g^) < cr) q2 (n)L q2n . (6.10) 

Finally we note that A(/i„, /2 ani0 ,o) — A(/x* , p an ,o,o)- Our next main aim is to estimate 
this quantity. 

By Proposition 3.1, G>*(z) = l/F^z), z G C + , where F^(z) := F^(Z(y/nz))/y/n. 
Here Z(z) G T is the solution of the equation (3.11) with /i = p* . 
Consider the functions 



S(z) := -{z + 5 n (z) := Z(yfcz)/y/K, 

Sni(z) a n + ^(z - a n + y/(z - a n ) 2 - 4j , 2 G C + . 

Note that 1/S(z) = G tlw (z), where /i w denotes the semicircle measure. Since G J 7 , we 
saw in Section 3 that there exists p n G M. such that l/S n (z) = Gp, n (z). In addition, it is 
easy to see, that 1/S nl (z) = G^^z). 
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In order to estimate A(/i* , // an ,o,o) we will apply the Stieltjes-Perron inversion formula 
to the measures \x* n and AV,,o,o- F° r this we need further estimates for | G M * (z) — G^ an (z) \ 
onC+. 

6.2. The functional equation for the function S n (z). Using (2.7) with /x = //*, we 
write, for z G C + , 

{ z(z)) = i + _ + _ y K 

= 1 + ^T + #T + ^T / (6-11) 



Z3(z) Z*(z)J R Z(z)-u 
The equation (3.11) with fJ, — /i* may be rewritten as 

G^{Z{z))[z{z) -z)=(n- 1)(1 - Z(z)G^(Z(z))), z G C+. (6.12) 

By (6.11) and the definition of S n (z), we represent (6.12) in the form 

1 , 1 r mWVon ^ n_1 1 fi , m 3 + ^n 2 (^)\ , R1 „v 

for z G C + , where 

r „ l(2 ) := / »VW , r-(t) := / |f + m ; - ros . (6.14) 



Z(y/nz) - u "~ J R Z(y/nz) - u 

By (6.1) and (6.4), we obtain from Lemma 3.4 for pL — pL* the bound 

\Z(Vnz)\ > y/(n-l)/8, z G C+. (6.15) 

The functions r n j(z), j = 1,2, are analytic on C + and with the help of the inequalities 
(6.3)-(6.6) and (6.15) admit the estimates, for z G C + , 



f \u\ 6 u*(du) 52/3* 53 /- 4 . Xr \ 

Vm{z)\ < / — L < -p. < — U + %2 (n)L q2n ) < 54L 3n , 

bw) < / j-rr + m 3 - m 3 < — ^- + 2^nr] q 2{n)L q2n 



< — p=r- + ^Vnrj q2 (n)L q2n . (6.16) 

We deduce from (6.13) the following relation 

Sl(z) - zS 2 n (z) + (1 + e nl (z))S n (z) + e n2 {z) = 0, z G C + , (6.17) 

where 

e n i(z) := -r„i(z) and e n2 (2:) := + r n3 (2) := a n + r n3 (2:) (6.18) 
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with 



r„ W := (l-i)^-~(l + r„ lW )- 
V nJ \ n n\ J 



m 3 



y/n n\ ) riy/n 

6.3. Estimates of e n i(z) and e n 2(z). By (6.16), we obtain 

\r n3 (z)\ < ^1 + 3r) q2 (n)L q2n + M(i + 54L 3n ) + ^, z G C+. (6.19) 
n n V / n 

Note that m 4 < /^(^n^ 4-92 ^ 2 . By (6.16), (6.19) and the last inequality, we have, for 
z e Di := {z e C + : < < 3, \^z\ < 4}, 

|e»iWI<54^<l (6.20) 
< 53 + 3 %2 (n)L g2n + 2kl+L3 " < 56 %2 (n)L, 2 „ + - (6.21) 

77/ 77/ 

and 

| £n2 (z)| < 56 %2 (n)L g2n + 2L 3n < (6.22) 

6.4. Roots of the functional equation for S n (z). For every fixed z G C + , consider 
the cubic equation 

P(z, w) := w 3 - zw 2 + (1 + e nl (z))w + e n2 (z) = 0. 

Denote roots of this equation by Wj = Wj(z), j = 1, 2, 3. 

We shall show that for z G D x the equation P(z,w) = has a root, say w 1 = Wi(z), 
such that 

wi = -a n + r n4 (z), (6.23) 

where the quantity r n4 (z) admits the following bound 

\r nA (z)\ < 10 2 r, where r := rj q2 (n)L q2n + L\ n . (6.24) 

In addition \wj + a n \ > 10 2 r, j = 2, 3. 
Indeed, introduce the polynomials 

Pi(w) := w 3 - zw 2 and P 2 (w) := (1 + e nl (z))w + e n2 {z) = (1 + e nl (z))(w + e n3 (z)), 

where e n3 (z) := e n2 (z)/(l + £ n i(z)). They admit the following estimates on the circle 
\w + a n \ = 10 2 r 

|PiH| < \w\ 2 (\w - z\) < 2{\w + a n \ 2 + a 2 n )(\w + a n \ + \z + a n \) 

< 2(10 4 r 2 + a 2 )(r + 11/2) < 12(10 4 r 2 + a 2 ) < 24 r (6.25) 

and 

|P 2 H| = \l + e nl (z)\\w + e n3 (z)\ > (1 - \e nl (z)\)\10 2 r - \e n3 (z) - a n \\. (6.26) 
Since, by (6.20), 1 - \e nl (z)\ > 9/10 and, by (6.20)-(6.22), 

\e n z{z) - a n \ < \r n3 (z)\ + 2je„i(^)| je r » 2 (^)| < 57r, 
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we see from (6.26) that | (xu) | > 36 r. This estimate and (6.25) gives us the inequality 
\Pi(w)\ < \P 2 (w)\ on the circle |w+a„| = 10 2 r and the desired result follows from Rouche's 
theorem. 

6.5. The remaining roots w 2 and w 3 of the equation (6.17). Now we shall inves- 
tigate the behavior of the roots w = W2(z) and w = w 3 (z). As seen in Subsection 6.4 
w 2 (z) 7^ Wi(z) and w 3 (z) 7^ Wi(z) for z G D\. We shall construct a set D 2 C D 1 where 
w 2 (z) 7^ w 3 (z), z G D 2 . Since P(z,w) = P 3 (z,w)(w — Wi), where 

P 3 (z, w) := w 2 - (z - wi)w + 1 + e nl (z) — Wi{z — Wi), 

we see that w 2 = w 3 for z G Di such that 

(z - Wl f - 4(1 + e nl (z) - Wl (z - Wl )) = 0. (6.27) 

We conclude from this relation that z = ±2-^/1 + e n i(z) + wf — W\. Therefore, as it is 
easy to see from (6.20) and (6.23), (6.24), the relation (6.27) does not hold for z G D 2 , 
where D 2 := {z G C : < ^sz < 3, \lftz — a n \ < 2 — hi} and h\ := c[ X ^r. 

Hence the roots Wi(z),w 2 (z) and w 3 (z) are distinct for z G D 2 . 

Now we see that the roots w 2 and w 3 have the form 



\(z-w l + (-iy- l ^Kz))., j = 2, 3. (6.28) 

where g(z) := (z — wi) 2 — 4 — Ae n i(z) + Aw\(z — w\) 7^ for z G D 2 . In this formula we 
choose the branch of the square satisfying yj g(i) G C + . 
Using (6.23), we rewrite (6.28) in the following way 

Wj := ^(z + a n + (-ly^y/iz + a n ) 2 - 4 - Aa n (z + a n ) + r n5 (2) j - ^r„ 4 (z) 

= a n + ^(z - a n + (-TT 1 y/(z - a n ) 2 - 4 - 4a 2 n + r n5 (z)^j - ^r n4 (z), (6.29) 
j = 1,2, where 

r n b(z) ■= -4e„i(z) + (2z + 6a n - 3r n4 (z))r n4 (z). 
From (6.20) and (6.24) it follows that the following estimate holds, for z G D 2 , 

\r n h{z)\ < 4|e„i(z)| + (2\z\ + Q\a n \ + S\r n4 (z)\)\r n4 (z)\ 

< 216— + (10 + 6L 3n + 300r)10 2 r < 1100 r. (6.30) 



Using (6.30) we obtain 



4a 2 - r n5 (z) 



1104 r i/ 6 1 , 

< — — < 1104 c^ 6 < — , zeD 2 . (6.31) 
hi 11) 



(2;-a„) 2 -4 

By power series expansion of (1 + z) 1 ^ 2 , \z\ < 1, we obtain, for z G D 2 

r n &(z) 



yj{z - a n ) 2 - 4 - 4a 2 + r n5 (z) = yj{z - a n ) 2 - 4 + 



v/(z-a n ) 2 -4' 
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where |r„6(z)| < 1004r. By this relation, we see that, for z G D 2 , 

Wj = a n + U(z - o„) + (-ly^V^-an) 2 -^ 

~ \r n ,{z) + ( " 1)J-1 rM , j = 2, 3. (6.32) 
2 2 v /( 2; _ an )2_4 

Let us show that S^z) = w$(z) for 2; G £>2- By (6.1), (6.23) and (6.24), we see 
that |wi(z)| < 1/6 for z G D 2 . Since S n (z) satisfies the equation (6.17) and, by (6.15), 
|>Sn(-2)| > 1/3 for all z G C + , we have S n (z) = w 2 (z) or S n (z) = w 3 (z) for z G D 2 . 

First assume that, for every z G D 2 , there exists r = r (-2 ) > such that S n (z) = 
Wj(z) for all z G D 2 fl { 1 2 — I < r o}, where j = 2 or j = 3. From this assumption it follows 
that S n (z) = Wj(z) for all z G -D2, where j = 2 or j = 3. Furthermore, it is not difficult 
to see that the roots w 2 (z) and w 3 (z) admit the estimates: ^(z)! < 4/3 for z G D 2 , and 
1^3(^)1 > 3/2 for z E D 2 and > 2. The analytic function S^z) G J 7 , by (3.6), satisfies 
the inequality SsS n (z) > ^sz. Hence, under the above assumption S n (z) = w 3 (z) for all 
z G D 2 . 

Now assume that the above assumption does not hold. Then there exists a point 
zq G D 2 such that, for any tq > 0, there exist points z' G D 2 and z" G D 2 in the 
disc \z — zq\ < r such that S n (z') = w 2 (z') and S n (z") = ws(z"). Let for definiteness 
S n (zo) = ""^(^o)- By this assumption, there exists a sequence {zk}^ =1 such that Zk — > z 
and S n (z k ) = w 3 (z k ). Therefore we have w 2 (z ) = lim 2fe _^ zo w 3 (zk). Using (6.32), rewrite 
this relation in the form 

a n + \ (z -On - V (zq - a n ) 2 - 4) - \ (r n4 (z ) + —=I^k^==) = 
2 V ) 2 V v /( Zo _ an )2_ 4 ; 

a n + J f^o ~a n + a/ (*o - a n ) 2 - 4) + J lim (r nA (z k ) r " 6 ^^ . (6.33) 

2 V / 2 z k -+zo V v /(^-a„) 2 -4/ 

From (6.33) we easily conclude with r as in (6.24) 

c' l lll2 ^-<W{z () -a n Y-A\ < 1004(c} /12 v^ + r), 

a contradiction for sufficiently small C\ > 0. Hence, the first assumption holds only and 
S n (z) = w 3 (z), z G D 2 . 

Denote by B x the set [—2 + /i x + a n , 2 — fox + a n ]. 

6.6. Estimate of the integral / |G> an (x + ie) — G M * (x + dx for < £ < 1. We 

Bi 

obtain an estimate of this integral, using the inequality 

/ \G^ anfifi (x + ie) -G^ n (x + ie)\dx < / \G^ anfifi (x + ie) - G^x + ie)\ dx 

J Bi J B\ 

+ / IG^ix + ie) -G^ n (x + ie)\dx. (6.34) 

JBi 
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Therefore we need to evaluate the functions 

G n an ,o,o ( z ) ~ G An (z) and Gjx n (z) - G„* (z) 

for z G D 2 . 

For z G D 2 , using the formula (6.29) with j = 3 for S n (z), we write 

1 1 S nl (z) - S n (z) 
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(6.35) 



G[i n {z) - G^ anfifi {z) 
r r ni (z) 4 



S n (z) S nl (z) S nl (z)S n (z) 
4a 2 n -r n5 (z) 



2S nl (z)S n (z)\ ^(z _ a n y -4+y/(z- a n ) 2 - 4 - 4a 2 + r n5 (z) 

By (6.31), we have, for z G -D 2 , 
\^/(z-a n ) 2 -A + y/(z - a n ) 2 - 4 - 4a 2 + r n5 (z) \ 



y (6.36) 



HV(*-On) 2 -4||l + v /l-(4a2-r n5 (z))/((z-a n ) 2 -4)| > | y/(z - a n ) 2 - 4|. (6.37) 

In addition, we see from (6.15) that (^(z)) > 1/3 for z G C + . The same estimate 
obviously holds for | ) S' n i(2;)|. 

Therefore we can conclude from (6.36) and (6.37) that 

f f 1 1 

/ G fin (x + ie)-G ltani0i0 (x + ie) 



dx 



S n (x + ie) S n i(x + is) 



dx 



<£ 




\r nA (x + ie)\ + 



4a 2 + \r n5 (x + ie)\ 
,/(x - a n + ie) 2 - 4 



dx 



(6.38) 



for < e < 1. 

From (6.24) it follows at once that 



/ |r n4 (x + i£)|dx<4-10 2 r, £ G (0, 1]. 



From (6.30) we conclude that, for the same e, 

, ^n + Mx + ie)] dx<1104r 
'b 1 W{x- a n + ie) 2 -4| 

It follows from (6.38)-(6.40) that 



dx 



< 4416 r. 



Bi 



G^x + ie) - G, Jx + ie 



Now we conclude from (6.11) that 
where 



Bi a/4 - (x - a n ) 2 

rfx<3-10 4 r, £G(0,1]. 



^(z) 



, -2 G C + , 



/ \ 1 r nl{z) 

r n 7\z) := ^77-^ + 



nSl(z) nS 2 (z) 



(6.39) 

(6.40) 

(6.41) 

(6.42) 
(6.43) 
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Since (^(z)! > 1/3 for z E C + , we see from (6.16) that 

K7W I < + z e D2 . (6 . 44) 

n 

Therefore, we deduce from (6.42) and (6.44) the upper bound 

/ \G^(x + ie)-G iin (x + ie)\dx< " 7V > dx < — — , ee(0,l]. (6.45) 
From (6.34), (6.41) and (6.45) we finally obtain 

/ \G M (x + te)-G^(x + te)\dx<A-10 i r, £6(0,1], (6.46) 

6.7. Application of the Stieltjes-Perron inversion formula. By (6.9), we have 
the relation 



Bi [2-^i+o„,2+a n ] [-2+a„,-2+Ai+a„] 

(6.47) 

From (6.46) and (6.47) we conclude, using the Stieltjes-Perron inversion formula, 

/<(#!) > 1 - (4 • 10 4 + q 1/4 r 1/2 ) r > 1 - (4 • 10 4 + c{ /4 ) r > 1 - (4 • 10 4 + 1) r. (6.48) 
Finally we deduce from (6.46)-(6.48) and the Stieltjes-Perron inversion formula that 

//a„,o,o) < cr = c (j] q2 (n)L q2n + L 2 n j. (6.49) 

6.8. Completion of the proof of Theorem 2.3. The statement of the theorem follows 
immediately from (6.7), (6.8), (6.10) and (6.49). □ 

Proof of Corollary 2.4. It is easy to see that the assertion of Corollary 2.4 follows from 
Theorem 2.3 and from the following simple formula, for x E R and n E N, 

I^a n ,0,o((-O0,x)) - (J, w ((-00,X)) = -^j=(x 2 - l)p w {x) +C#(^=!) ' . 



Proof of Corollary 2.5. The assertion of Corollary 2.5 follows immediately from (2.20) 
and Proposition 3.6. □ 

7. Edgeworth Expansion in Free CLT (the case (3 q < oo, q > 4) 

In this section we prove Theorem 2.6 and Corollary 2.7. The proof of the theorem is 
similar to the proof of Theorem 2.3 but with some essential technical differences. Therefore 
we describe in detail those arguments which differ from the proof of Theorem 2.3 and omit 
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arguments which directly repeat the arguments of Section 6. We preserve all notations of 
Section 6. Denote as well 

S n2 (z) := a n + + b n ){z - an) + sj (1 - b n )\z - a n f - 4(l - d n )), z £ C + , 

where a n , b n and d„ are defined in Section 2. The function S n 2(z) £ J 7 and l/S n 2(z) = 
bn dn (z), where fj, an ,b n ,d n is the free Meixner measure with the parameters a n , 6 n and 
d n , see (2.14). 

Proof of Theorem 2.6. First we proceed to study 



7.1. The passage to measures with bounded supports. Let n £ A/". Let e n £ 
(0, 10" 1 / 2 ] be a point at which the infimum of the function g qn z{s) in (2.16) is attained. 
This means that 

vAn) :=el^ + Pq3 ^ R £nVE) e^. (7.1) 

Using this parameter e n , we define free random variables X, X 1: X 2 , . . . and X*, X*, X%, . . . 
in the same way as in Section 6. We define probability measures /i n , //*,//* in the same 
way as well. 

Without loss of generality we assume that 

r) q3 (n)L q3n + L 4n < c 2 , (7.2) 

where c 2 > is a sufficiently small absolute constant. From (7.2) it follows that n is 
sufficiently large n > c^m^ > c^ 1 - Here and in the sequel we use Lyapunov's inequality 
1 = m^ 2 < (3^ 3 < m^ 4 . 

Now we repeat the arguments of Subsection 6.1. 

Using (7.2) we note that 

14,1 < e^-Vn-te-Wpnfaeny/E) < ^= Vq ,{n)L qsn (7.3) 

and 

< -k- ~ 1 < 2(p 2 (/i, e«v^) + K) < 3r) q3 (n)L q3n . (7.4) 
By (7.3) and (7.4), we see that (6.4) holds and the support of ji* is contained in [— \y/n, \\fn\ 
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Recalling (6.5) and (6.6), we easily deduce, by (7.2)-(7.4), that 

\m* 3 ~ m 3| < 2y/nr] q3 (n)L q3n , 
\m* A — m 4 | 

< C~ 4 \m A - m 4 1 + (C~ 4 - l)m 4 + C- 4 (4|A n ||m 3 | + 6A'm 2 + 4|A n | 3 |m 1 | + A 4 



< C n 4 \m 4 - m 4 \ + 5m 4 r) q3 (n)L q3n + 5 ^ + 1 r) q3 (n)L q3n 

< C^e-^n-^-^p^e^) + 5m 4 (l + -) Vq3 (n)L q3n 

<2nr] q3 (n)L q3n (7.5) 

and 

/3 5 * < C" 5 & + 6L 4nVq3 (n)L q3n . (7.6) 
By the triangle inequality we have 

A(//„, /t„) < A(/i„, /2 n ) + A(/i„, R n ) + A(/c n , k„), (7.7) 

where, for x G R, 

K n ((-oo,x)) := Pa n ,b n An{{-oo,x)) + -^,((-00, x)) 

■= Ha n ,b n ,d n ((-oo, (x - y/nA n )/C n )) + -<r n ((-oo, (x - y/nAn)/C n )). 

Th 

Note that A(/i n , it n ) = A(/x* , « n ). 

First we establish with the help of Proposition 3.5 

A(Hn,p> n ) < nA(fi,p,) < n/j,({\u\ > e n y/n}) < e~ qs n~ {q3 ~ 2)/2 p q3 (ii : e n y/n) < r] q3 (n)L q3n . 

(7.8) 

We saw in Section 4 that, for n > 3m 4 , p an ,b n ,d„ is an absolutely continuous measure 
with support on the set B 2 :— [a n — 2/e n , a n + 2/e n ] and density of the form 

_ y/4(l - jg - (1 - b n ) 2 (x - a n ) 2 
P» aMn W - + fln(1 _ fe> + 1 _ 4) ' U ' 9j 

This density does not exceed 1 on the set B 2 and is equal outside of this set. 

The signed measure q n has density p fn , see (2.25), which does not exceed 1 by modulus 
on the set B 2 and is equal to zero outside of B 2 . 

Therefore, in view of (7.2), a simple calculation shows that 

A I A 1 1 / 1 JnA n \ 

A(k„, K n ) < A(/x an)6 „ )dn , p an ,b n ,dj + -A(^ n , < c — - 1 + — — 

< ce-^-^n-^-^/^gs^, e«M < c Vq3 (n)L q3n . (7.10) 
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Our next aim is to estimate the quantity A(//*, K n ). In order to estimate A(//*,/c n ) 
we need to apply the inversion formula to fi* n and K n . We shall now derive the necessary 
estimates for \G^(z) — G Kn (z)\ on C + . 

7.2. The functional equation for S n (z). Using (2.7) with fi = fi*, we write, for z G C + , 

, WG , (ZW) . x + _L_ + + + ^/^. an, 

By (7.11) and the definition of S n (z), the equation (3.11) with fx = fx* may be rewritten 



as 



^ 1 m 3 . rn* 4 + ( nl (z) \ 

\ + nS*(z) + nWS*(z) + n*S*(z) ) [n[) ' 

_ n-l 1 / ^ m* | m* + Cni(^) \ , ? ^ 

for where C„i(*) := / R f^} u - 

We deduce from (7.12) the following relation, for z G C + , 

S*(z) - zS*(z) + S*(z) + C -^Sl(z) + C -^S n (z) - = 0, (7.13) 

\f 71/71/ 71/ 

where C, n i{z) := ™\ - z/y/n, ( n3 (z)(z) := m* + ( nl (z) - zm* 3 /y/n and C,m{z){z) : = m* + 
(ni(z). Note that the functions ( n j(z),j = 1,2,3,4, are analytic on C + . 

7.3. Estimates of the functions ( n j(z), j = 1,2,3 on the set D x . From (7.2) and 
(6.4) we deduce, using Lemma 3.4, that (6.15) holds. Therefore, in view of (7.4) and 
(7.6), we arrive at the estimate 

U7(rM | || . ^ < 53 A + 312 ^vML, m 

M<§>/n \\ Z W nz )\ ~ \ U \\ V n V n V n 

<53 Pq3{ n } + 312-^J7,3(n)^n<54n77, 3 (n)L wn , z G C + . (7.14) 

For z G Di, by (7.5) and (7.14), we get the bounds 
K " 2(:)l <2(L 3n + Vq3 (n)L q3n ), 





u\ 


5 fi*{du) 


\\Z(yfiz)\-\ 


u\\ 



n 



£ m ; + M.-)i + ski £ 2Ltn + 56th3(n)Kn 



n n n 



\Cn4z)\ .ml+\( nl (z)\ 

< < L 4n + 56rj q3 (n)L qsn . (7.15) 



n n 
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7.4. The roots of the functional equation (7.13) for S n (z). For every fixed z G C + 
consider the equation 

r\l \ 5 4,3, Cn2\z) 2 , Cn2,i.z) ( ni (z)Z 

Q(z, w) :— w — zw + w H H w ^ — = 0. (7.16) 

y/n n n z 

Denote the roots of the equation (7.16) by Wj = Wj(z), j — 1, . . . , 5. Let us show that 
for every fixed z G D\ the equation Q(z,w) = has three roots, say Wj = Wj(z), j = 
1, 2, 3, such that 

\ Wj \ < r' := 15(L 4n + %3 W^ 3n ) 1/2 , j = h 2, 3, (7.17) 

two roots, say Wj, j = 4,5, such that \wj\ > r' for j = 4,5. Recalling (7.2) we see that 
the bound r' < holds. 
Consider the polynomials 



n ( \ 5 4 , Qi\ Z ) 2 , Cn?,{z) Cni{z)z 3 

Qi{z, w) :— w — zw H )=-w A w — and Q 2 (w) :— w 

Jn n n 2 



The following estimates hold on the circle = r' for z G Di. |u>| 5 = (r') 2 \w\ 3 < 10 _4 |w| 3 , 
and \zw A \ = |z|r'|w| 3 < ^j|w| 3 - Since nr' > lhy/ml > 15m 2 = 15 and, by Proposition 3.6, 
r' > 15L 3n , we have as well, using (7.15), 

lCn2(*)| ,2. = \U{Z)\ 1 3 < |uf 



y/n y/n r' 5 

ICnsWL , \U(Z)\ 1 3 A\ W \ 

\w\ = — -to < 



n n (r') 2 15 

ICn4(^L, _ |Cn4(*)| \z\ , ,3 < 4|*| 3 < 4[ W | 



n 2 n 2 (r') 3 15w' 45 

We see from the last five inequalities that \Qi(z, w)\ < ^\Q 2 (w)\ on the circle \w\ = r' . 
Therefore, by Rouche's theorem, we obtain that the polynomial Qi(z,w) + Q 2 (w) has 
only three roots which are less than r' in modulus, as claimed. 
Represent Q(z,w) in the form 

Q(z, w) = (w 2 + siw + s 2 )(w 3 + giw 2 + g 2 w + g 3 ), 

where w 3 + giw 2 + g 2 w + g$ = (w — wi)(w — w 2 )(w — W3). From this formula we derive 
the relations 

Cn2{z) 



si+gi = -z, s 2 + si^i + g 2 = 1, s 2 g x + sig 2 + #3 



S 2 #2 + Si5- 3 = , S 2 g 3 = 5 • (7-18) 

n n z 

By Vieta's formulae and (7.17), note that 

|<7i| < 3r', \g 2 \ < 3(r') 2 , |<? 3 | < (r') 3 . (7.19) 
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Now we obtain from (7.18) and (7.19) the following bounds, for z G D\, 

|si| < 5 + 3r, |1 - s 2 \ < 3r(4r + 5) < 16r < ^. 
Then we conclude from (7.5), (7.15), (7.18)-(7.20) that, for the same z, 
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(7.20) 



92 ~ 



n 



< 



92 



Cn3{z) 



n 



., ; .<W w + fe_HM + (r , f 

n 6 ^ \S2\ \S2\ n 



< n{r'Y + S(r'Y + (r'Y = 20(r') d . 



(7.21) 



Denote a* := m* 3 /y/n, p* n := L* 4n -l/n := (m* 4 -l)/n and p n := (m 4 - l)/n . By (7.5), 
it is easy to see that 

\a n - a* | = y=-^- < 2rj q3 {n)L q:i , n and \p n - p* n \ = < 2rj q3 {n)L q . }n . (7.22) 

From the first three relations in (7.18) it follows that 

9i + zgf = a n + p n z + ( n5 (z), (7.23) 

where 



Cns(z) := ^L^l + L 2 - kli^TW -gf + 2 gi g 2 -g 3 - (a n - a* n ) - (p n - p* n )z. 
n V n / 

By (7.14), (7.19), (7.21) and (7.22), we get the following estimate, for z e D u 



\U{z)\<\^^ + 



n 



92 



n 



\A + \9i\ +2|ffi^ 2 | + M + K~<| + |p„-p;| 



< 274 Vq3 (n)L qsn + 146(r') 3 < 8 • 10 5 (r] q3 {n)L qzn + L%^. 
Rewrite (7.23) in the form 

+ a n z) = a n + p n z + (a n + p n z) ( 1 - - l) + a n #iz + 



(7.24) 



Taking into account (7.19), (7.24) and Proposition 3.6 this relation leads us to the bound, 
for z E D\, 



\9i ~a n - {pn - a 2 n )z\ 



< \ a l z<1 \ 



\Q"nPnZ , 



I 2 2 1 
|On0i2 I 



+ 



+ 



|l + a„z| \l + a n z\ |l + a n 2||l + ^i2;| |1 + a n z\\l + g x z\ |1 + a n z\\l + g x z\ 
< 50L 3 3n + 50L 3 „L 4ri + 500L3„(r') 2 + 150L 4n r' + 16 • 10 5 (r] g3 (n)L q3n + L^ 2 ) 

<18-10 5 (^ 3 (n)L g3ri + L^ 2 ). (7.25) 
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To find the roots W4 and W5, we need to solve the equation w 2 + siw + s 2 = 0. Using 
(7.18), we have, for j = 4, 5, 

= \ (z + gi + y/(z + 9i) 2 - + + 9i)9i -9*)) 

= I (z + 9l + (-iy J ( Z - 9l )* - 4 - 4(g* - frj) = ^U(z) + a n + 

+ X - ((l + b n ) (z - a n ) + (-1)' ^(1 - b n )\z - a n ) 2 - 4(1 - d n ) + C„ 7 (z)) , (7.26) 

where 

Cne(z) := gi-a n - b n (z - a n ), 

( n7 (z) := -3£ 6 (z) - 2U{z){4a n + (1 + 3b n )(z - a n )) 

+ 4(# 2 - L 4n ) - Ab n (z - a n )(2a n + b n (z - a n )). (7.27) 

We choose the branch of the analytic square root according to the condition ^sw 4 (i) > 0. 
Note that the roots w±{z) and w 5 (z) are continuous functions in D\. 

7.5. Estimates of the functions ( n6 (z) and ( n7 (z) on the set Di. We obtain, by 
(7.25), 

\U(z)\ < 18-10 5 (^ 3 (n)L g3n + L^ 2 ) +|a n 6 n | < (18 • 10 5 + 1) (vM^n + L 3 ^ . (7.28) 
By (7.14), (7.21) and (7.22), we have 

\g 2 - L 4n \ < 56r] g3 (n)L q3n + 20(r') 3 < 10 5 (r] q3 (n)L q3n + (7.29) 
Then, using (7.2), (7.28) and (7.29), we easily deduce from (7.27) 

\Cn 7 (z)\ <3|Cn 6 W| 2 + 2|C„6(*)|(4|a n | + (1 + Sb n )(\z\ + \a n \)) + A\g 2 - L 4n | 

+ 4b n {\z\ + \a n \)(2\a n \ +b n {\z\ + \a n \)) < 3 • 10 7 (r] q3 {n)L q3n + L^ 2 ) . (7.30) 

7.6. The roots W4 and w$. We saw in Subsection 7.4 that w±(z) 7^ Wj(z), z G Di, for 
j = 1,2,3. Returning to (7.26), it follows from (7.30) that Wi{z) 7^ w${z) for z G D 3 , 
where 

D 3 := [z G C : < %z < 3, \Uz - a n \ < h 2 \ 

where e n := (1 - b n )/y/l - d n and /i 2 := C2 1/6 (j) q3 {n)L q3n + L^ 2 j . 

Since the constant c 2 > is sufficiently small, we have, by (7.30), for z G D 3 , 

|Cn7(«)|/|((l - b n ) 2 (z - a n f - 4(1 - d n )\ < 4 • 10 7 cf < lO" 2 . (7.31) 
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Therefore, using power series expansion of (1 + z) 1 ^ 2 , \z\ < 1, we obtain, for the same z, 
- b n ) 2 (z - a n ) 2 - 4(1 - d n ) + Cm(z) = y/(l - b n ) 2 (z - a n ) 2 - 4(1 - d n ) 

v /(l-6„) 2 (^-a n ) 2 -4(l-rf„)' 
where |Cns(z)| < 4 • 10 7 (j] q3 (n)L q3n + L^f^. By this relation, we see that 

Wj 



(z) = a n + -[(l + b n ) (z - a n ) + (-iy \j (1 - b n ) (z - a n ) 2 - 4(1 - d n 

+ -U(z) + ; Cw8( ^ j = 4, 5, (7.32) 

2' 61 ) 2 ^(l-^- ^ -4(1-0 

for 2 G £> 3 . 

Let us show that S n (z) = w±{z) for z G -D3. By (7.2) and (7.17), we see that |iu,-(;z)| < 
1/6 for z G D 3 . Since, in view of (6.15), (^(z)! > 1/3 for all z G C + , we have S n (z) = 
Wi{z) or S n (z) = w$(z) for z G D 3 . 

Assume that, for every z G D 3 , there exists r = ro(zo) > such that S n (z) = Wj(z) 
for all z G D3 H {\z — z \ < r}, where j = 4 or j = 5. From this assumption it 
follows that S n (z) = Wj(z) for all z G D 3 , where j = 4 or j = 5. By (3.6), we have 
\S n (2i)\ > I. In addition it follows from (7.17) and (7.26), (7.28), (7.30) that \w A (2i)\ > 1 
and \wj(2i)\ < 1, j — 1, 2, 3, 5. Hence in this case S n (z) = w±{z) for z G -D3. 

If the above assumption is not true, there exists a point z G D 3 such that, for any 
r > 0, there exist points z' G D 3 and 2" G -D3 from the disc \z — z \ < r such that 
S n (z') = w±{z') and S n {z") = w 5 (z"). Let for definiteness S n (z ) = w 5 (z ). By assumption 
there exists a sequence {z k }^ =1 such that £ fc — > z and S n (zk) = w±(z k ). Therefore we 
have ws(zo) = ]im Zk ->. Zo w±{zk)- Rewrite this relation, using (7.32), 



a n + ^((l + b n )(z - On) - b n ) 2 (z - a n ) 2 - 4(1 - d f 

+ 1 t „ \ Cns(^o) 



V (z)- Us^oj \ = 

^ v/(l-M 2 (^o-an) 2 -4(l-rf„)^ 

On + ^ ((l + K) (zo - a n ) + \j (l- b n ) 2 (z - a n ) 2 - 4(l - d n ) j 
+ i lim f U(^) + Ui ~ Zk) =) ■ (7.33) 

From (7.33) we easily conclude 

c 2 1/12 \/ Vq3(n)L q3n + L 3 4n 2 < - b n ) 2 \z - a n ) 2 - 4(1 - d n ) | 

< 4 • 10 7 y / ^ 3 (^)^g 3 r t + ^4n 2 + \J Vtfin) L q3n + , 
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which leads to a contradiction for sufficiently small C2 > 0. Hence our assumption holds 
and S n (z) = w±{z), z G D 3 . 

Denote by B 3 the set — - + h 2 + a n , - — h 2 + a n \. 

7.7. Estimate of the integral J Bs \G K (x + ie) - G^ bndn (x + ie) - ^{G„ aMn {x + 
ie)) 3 \ dx for < e < 1. We obtain an estimate of this integral, using the inequality 

/ \G^(x + is) - G^ anbndn {x + ie) - Ug^^x + ie)f\ dx 



< / \G„ aMn (x + ie) -G iXn (x + ie)\dx+ 

'B 3 



[ \G fln (x + ie)-G K (x + ie) - Ug^^x + ie)f\dx. 



(7.34) 



Therefore we need to evaluate the functions G^ an bn dn (z) — G^ n (z) and G^ n (z) — G M » (z) — 
~(Ga „ , {z)f for z G I),. 

For z G D 3 , using the formula (7.26) with j = 4 for S n (z), we write 



S n2 (z)S n (z) (— ^ - 1 ) = S n2 (z) - S n (z) = -\Cn&{z) 
\o n {z) o n2 {z)/ Z 

1 Cn7(z) 



2 - b n f{z - a n f - 4(1 - dn) + y/(l - b n )\z - a n f - 4(1 - d n ) + ( n7 (z) ' 

(7.35) 

Using (7.31) we get, for z G D 3 , 



y/(l - b n )\z - a n f - 4(1 - d n ) + y/(l - b n f{z - a n y - 4(1 - d n ) + Cm{z) 



^{l-b n y{z-a n f-A{l-d n ) 1 + y/1 + ( n7 (z)/((l - b n f{z - a n y - 4(1 - d n )) 
>W{l-b n f{z-a n f-A{l-d n )\. (7.36) 

It is easy to see that the bound |<S' n (2;)| > 1/3, z G C + , holds for S n2 (z) as well. Therefore 
in the same way as in the proof of (6.39) and (6.40) we conclude from (7.28), (7.30) and 
(7.35) that 



/ Gjx n (x + ie) - G^ aMn (x + ie) dx = 



1 



< 



B 2 



\( n6 (x + ie) \ + 



S n (x + ie) S n2 (x + ie) 
\( n7 {x + ie)\ 



\y/(l-b n )*(x-a n + ie)*-4(l-d n ) 



dx 
dx 



< c(r] q3 (n)L q3n + Ll^. 



(7.37) 
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Now we deduce from (6.11) with the probability measure //*, involving e n introduced 
in (7.1), the relation 

1 r nl (z) 



G^-G^z) 



-, ze C + , 



(7.38) 



nSl(z) nS*(z) 

where the function r n i(z) is defined in (6.14). Since (6.20) holds for r nl (z), we see that 



rm(z)\ <145g ^3n_ 



Since, for the same z, 

1 1 1 

< 2 



n\Sl{z 
1 



z G As- 



n 



S 3 n (z) Sl 2 {z) ~ S n {z) S n2 (z) \\S n (z)\* \S n2 (z)\ 



< 36 



S n (z) S n2 (z 



(7.39) 



(7.40) 



we obtain, using (7.37)-(7.40) and Proposition 3.6, that, for < e < 1, 



JB 3 






f 




B 3 


1 




+ - 


J 


n 






B 3 



dx 



dx 



(G» an , bn , dn (x + is)? ~ (G^x + ie)f\ dx < c {n q ,{n)L q , n + L% 2 ). (7.41) 



From (7.34), (7.37) and (7.41) we finally get, for < e < 1, 



\G u *(x + ie) - G 



b 3 



» aMn + is) - -(G^ n 6nidn (x + is)) | dx 



< c (r} q3 (n)L, 



■+- r 3/2 } 



(7.42) 



7.8. Application of the Stieltjes-Perron inversion formula. Using (7.9), we have 
the relation 

/ P, an , bn ,A*)dx>l-hl /2 . (7.43) 
Jb 3 

where hn dn (x) denotes the density bn dn (x) of the measure fi an ,b n ,d n It is not difficult 
to verify that 

q n i(dx) fP<; nl (x)dx 
Jm. 



z — X 



z — X 



-, ze C+, 



where 



X 



3((1 + b n )x + (1 - b n )a n f + (1 - b n ) 2 {x - a n ) 2 - 4(1 - d n ) 
(b n x 2 + (1 - b n )a n x + 1 - d n y 
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for x G B 2 and Pq nl {x) = for x B 2 . Therefore we easily deduce the obvious upper 
bounds 

/ Pfci(aj) dx < ch 2 /2 > / \P^i( x )\ dx < ch l /2 and A^nl^n) < c (\a n \ + L 4n ) . 
Jb 3 Jr\b 3 v 7 

(7.44) 

From (7.42)-(7.44) and the Stieltjes-Perron inversion formula we conclude that 

<(£ 3 ) > 1 " c (ifcsML^ + Lj?) . (7.45) 

We finally conclude from (7.42), (7.44), (7.45) and the Stieltjes-Perron inversion formula 
that 

A(/<, n n ) < c (j] q3 (n)L q3n + L 3 ^ . (7.46) 

7.9. Completion of the proof of Theorem 2.6. The statement of the theorem follows 
immediately from (7.7), (7.8), (7.10), (7.46) and Proposition 3.6. □ 

Proof of Corollary 2.7. Recalling the definition of the density p^ an bn dn {x) of the mea- 
sure Ha n ,b n ,d n for n > C2 1 m 4 , we see that 

P» an , bn , dn ( X + a n) 

= ^-y/(4(l-d n )-(l-b n )W) + (l + d n -b n - a n x - (b n - a 2 n )(x 2 - 1)) 
+ c6{L in + a 2 n ) 3/2 , xeM. 
In addition we have, for ieK, 

^ f ^(A{l-d n )-(l-b n ) 2 u 2 ) + du 

= (1 - d n + b n )fi w {{-oo,x)) + {-d n - a/(4 - x 2 ) + + c 0I?J? 

and, for x G (— l n ,l n ), where l n := max{2,2/e„}, 



| v/(4(l - 4) " (1 " b n ) 2 X 2 ) + " V 7 ^^) 



+ 1 

cZ/4 n 



v/(4(l - 4) - (1 - &„) 2 * 2 )+ + V(4 - * 2 ) + ' 
Using these formulae and the following obvious relations 

J xVA-x 2 dx = -i(4-x 2 ) 3/2 and ^ (x 2 - 1)VZ - x 2 dx = ~x(A - x 2 ) 3/2 , 

we obtain from (2.26), using some simple calculations, the representation (2.27). □ 
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8. Proof of Theorem 2.10 

In this section we keep the notations of Sections 6, 7. By definition of the random 
variable X* (see Subsection 7.1 in the case g 3 = 5) we note that \x* = pi for n > 10L 2 . 
Hence for these n we have //* = /i n as well. By Lemma 3.3, S n (z) is continuous up to the 
real axis and, by Lemma 3.4, l-S^z)! > 1.03/3 for z G C + U R and for n > c(/j,). Since 

c»(») = / , . f"\ r , 

J R SJz) -u Jn 



where supp(fi) C [-L, L], we conclude that G fin (z) is a continuous up to the real axis. By 
Lemmas 3.3, 3.4 and (3.4), fi n and ji n are absolutely continuous measures with continuous 
densities p n (x) and p n (x), respectively, defined by 

Pn{x) = -lim-5 — — and p n (x) = - lim -SG„ n (x + is). 

eiO 71 S n (X + te) eiO 7T 

In addition, p n (x) < 1 and p n (x) < 2 for all x G R and n > Ci(/x). We will use the 
estimate, for iGl, 

\p n {x)-p^ anibntdn {x)-^p, nl {x)\ < \Pn{x)-p^ an:bn:dn X)\^p n {x)-p n {x)-^ nl {x)\. (8.1) 

Repeating the arguments of Subsection 7.7, we easily obtain the upper bounds, for x G B 3 
and n > c(/i), 

T -i- r 3 / 2 

Iz5(a0-D I < c " 4n (8 2) 

a/4(1 - rf„) - (1 - b n y{x - a n y 

and 

\ Pn (x) -p n (x) - l -p, nl {x)\ < + j=^f7 ( 8 - 3 ) 

n V n nv /4(l - d n ) - (1 - 6 n ) 2 (x - a„) 2 / 

On the other hand, by the formula (7.9), we have 
P» an>bn , dn ( x ) = (i + ^d n -a 2 n -a n (x-a n )-(b n -a 2 J(x-a n ) 2 ^p w (e n (x-a n ))+c9Ll n 2 (8.4) 
for x G B 2 , where e n := (1 — b n )/y/l — d n . Returning to Subsection 7.8 we note that 

PwW = i( x - a nf - l)p w (e n (x - On)) + c9(\a n \ + L 4n ) (8.5) 

for x G B 2 . 

Applying (8.2)-(8.5) to (8.1) we arrive at the relation 

3/2 



where 



p n (x) = v n (x - a n ) + c9 Ls " ± xeB 3 , (8.6) 

a/4(1 - dn) - (1 - 6„) 2 (a; - a n ) 2 

«n(aj) : = (i + 7^ n -a 2 n -^- a n x - (b n -a 2 n - ^jx 2 ^jp w (e n x), i£l 
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The statement (2.31) of the theorem easily follows from (8.6). Thus, Theorem 2.10 is 
completely proved. 

9. Asymptotic expansion of / r \p n (x) —p w {x)\ dx as n ->■ oo 

In this section we shall prove Corollary 2.11. Note that for free bounded random 
variables the inequality (7.45) has the form /j, n (B s ) > 1 — c(L 5n + L% 2 ). Hence, we have 

I \p n (x) - p w (x)\dx = / \p n (x) - p w (x)\dx + 9(c\a n \ 3/2 + ^jjl). (9.1) 
Jr Jb 3 n 

First let us assume that m 3 ^ 0. Using (8.6), we easily conclude that 

\p n (x) - p w (x)\ dx = / |(1 — a n x)p w (x) — p w (x + a n )\ dx + 9[c\a n \ 3 ^ 2 H ) 



B 3 J [-2,2] 



2tt 



f \x\^z£idx + 9{c\a n f 2 + C M) = 2 ±A + d (c\a n f 2 + ^). 
J [-2,2] V4-x 2 n 7T n 



Applying this relation to (9.1) we arrive at the expansion (2.33). 
Let m 3 = 0. By (8.6), we easily get 



\p n (x) -p w (x)\dx 

B 3 

2. 

p w {x) 



I Pw{x) -p w (— ) + (]:d n - - - (b n - 

7f-2,2l Ve "^ V2 n V n/e n 



e„ n 3 / 2 



1 



2tT 7[-2,2] V4 - x 2 

| g cfc) _ [m 4 - 2[ y |2 - 4x 2 + x 4 | ^ + g c(/i) _ 2|m 4 - 2| | g cfc) 



n 3 / 2 27m 7[-2,2] V4 — x 2 ^ 3 ^ 2 7m ^ 3//2 

Thus, (2.34) is proved. 

10. Asymptotic expansion of the free entropy and the free Fisher 

information 

In this section we prove Corollaries 2.12 and 2.13. First we find an asymptotic expan- 
sion of the logarithmic energy E([i n ) of the measure \i n . Recall that (see [25]) 

-E(n n ) = / log \x - y\ fi n {dx)n n (dy) = h{fi n ) + h{^n) ■= 
J Jr 2 

= // log \x -y\fi n {dx)fi n {dy) + / / log \x - y\ fi n (dx)fi n (dy). 

J JB 3 xB 3 J JR 2 \(B 3 xB 3 ) 

(10.1) 
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Using (2.32) and the estimates p n (x) < 2, \x\ < 2 + L/y/n and p n (x) = 0, \x\ > 2 + L/y/n, 
we conclude that 

\h(l^ n )\<4= Pn{x) / \\og\x - y\\p n (y) dy dx < c(n) / p n (x)dx 
Jr\b 3 Jr Jr\b 3 

< c{fi)n- 3 / 2 . (10.2) 
By (8.6), we have 




h{Hn)= I I log|x - y\v n (x - a n )v n (y - a n )dxdy + c(/i)9(L 5n + L^ 2 ) 

J JB 3 xB 3 

= 11 log \x -y\v n (x)v n (y)dxdy + c(fi)6n- 3/2 . (10.3) 
Recalling the definition of v n (x) we see that 

log | a: - y\ v n (x)v n (y) dxdy = h{v n ) + h{v n ) + h{v n ) + h(v n ) + r n 

:= (l + \d n - a 2 n - -) I I log \x -y\p w (e n x)p w (e n y) dxdy 
v ^ n/ J jr2 

- 2^1 + \d n - a 2 n - -)a n I I x\og\x - y\p w {e n x)p w (e n y) dxdy 

v 1 n/ J JR2 

+ a 2 n xy log \x -y\p w (e n x)p w (e n y) dxdy 

J JR.* 

- l{b n - a 2 n - -) I I x 2 \og\x - y\p w {e n x)p w (e n y) dxdy + c([i)9n~ 3/2 . (10.4) 



In view of E(n w ) = 1/4 and e n 2 = 1 - d n + 2b n + c9L\ n , log e n = 4f - b n + c9L\ n , note 
that 



h(v n ) = (l + ^d n - a 2 n - e « 2 ( / Jog\x -y\p w (x)p w {y) dxdy- log e n ^j 



= -E(^) + ^ + c(/i)^ 2 . (10.5) 

Since the function J R log|:r — y\p w (y) dy is even, we see that hivn) = 0. In order to 
calculate h(v n ) we easily deduce that 

J up w {u) log \x — u\ du = — x + x 3 /6, x & [—2,2]. 

Therefore we obtain 

/ 2 2 
xp w (x)(x - x 3 /6) dx + c(fi)9n- 2 = --a 2 + c(/i)#n~ 2 . (10.6) 
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Using the following well-known formula (see [25], p. 197) 



x 2 1 

p w (u) log \x — u\ du — — , x G [—2, 2], 



we deduce 



J J ^x 2 log \x -y\p w (x)p w (y)dxdy = J x 2 - ^jp w (x) dx = 0. 

Therefore 

h(v n ) = c(fi)6n- 2 . (10.7) 
By (10.3)-(10.7), we arrive at the formula 

h{ix n ) = -E{fi w ) - ^a 2 n + c{fi)6n- 3/2 . 

Finally, in view of (10.1)-(10.2), we get 

-E(fi n ) = -E(n w ) - l -a 2 n + c(fi)9n-^ 2 

and we obtain the assertion of Corollary 2.12. 

Now let us prove Corollary 2.13. We shall show that the free Fisher information of the 
measure \i n has the form (2.35). Denote 

4tt 2 f 4rr 2 f 

= $iW + $ 2 (/v) := -r- / p n (xfdx + — Pn{xfdx. (10.8) 

-J Jb 3 -j Jr\b 3 

As before we see that, by (2.32), 

$2 W < c(/i) / Pn (x) dx < c(/i)n- 3 / 2 . (10.9) 

JR\B 3 

On the other hand, by (8.6), we have 

®M = / u„(a:) 3 da; + c(Ai)0n- 3 / 2 . (10.10) 
It is easy to see that the integral on the right hand-side of (10.10) is equal to 

(l + ?>0-d n - a 2 n - -Yje' 1 / p w (xfdx - ?>(b n - a 2 n - -W 3 / x 2 p w (xf dx 

+ 3 «n e n 3 / ^V(^) 3 ^ + c(fi)0L%? = — (1 + a 2 ) + c(n)9n- 3/2 . 

Therefore we finally have by (10.8)-(10.10) 

$(/i„) = 1 + a 2 + c(fi)9n- 3/2 = $(/i w ) + a 2 + c(/i)#n~ 3/2 . 
Thus, Corollary 2.13 is proved. 
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11. Appendix 1. Proof of Theorem 2.1 
In this Appendix we keep the notations of Section 6. 
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11.1. Passage to measures with bounded supports. Let n G Af. Let e n G (0, 1CT 1 / 2 ] 
be a point at which the infimum of the function g qn i(s) from (2.16) is attained. This means 
that 

Pqi 

Using this parameter e n , we define free random variables X, Xi,X2, ■ ■ ■ and X*, X*, X%, . . . 
in the same way as in Section 5. We define probability measures p n , p w = po,o,o, /•**, m 
the same way as well. 

Without loss of generality we assume that 

r} ql (n)L qin + 1/n < c 3 , (11.1) 

where C3 > is a sufficiently small absolute constant. 
Using (11.1) we note that 

14,1 < e-^n-^-^p^e^) < -^=r) ql {n)Kn (H-2) 

and 

< ^- - 1 < 2{p 2 { l i,e n y/n) + A 2 ) < 3 Vql (n)L qin . (11.3) 

By (1 1.1)— (11.3), we obtain that (6.4) holds and the support of p* is contained in 
[—\\/n, \^/n\- By (11.1)— (11.3), we easily deduce as well that 

^<C- 3 ^ + ^= Vql (n)L qin . (11.4) 

By the triangle inequality, we have 

A(p n , p w ) < A(p n , p n ) + A(p n , p w ) + A(p w , p w ). (11.5) 

Furthermore, we have the following inequalities 

A(p n ,p n ) < e- qi n- {qi - 2)/2 p qi (ii,e n Vn) < r] ql (n)L qin , 

A(p w ,p w ) < ce-^-^n-^- 2 ^ 2 p qi (p,e n V^) < cr) ql (n)L qin . (11.6) 

Our next aim is to estimate A(p n ,p w ) = A(p* n ,p w ). 

As in Section 6, let Z(z) G T be the solution of the equation (3.11) with p = p*. 
Denote S n (z) := Z(y/nz)/ 'y/n. 
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11.2. The functional equation for S n (z). Using the formula 

z '*'»» = 1 ^ + ^/,S (1L7) 

and the equation (3.11) with /i — /i* we arrive at the following functional equation for 

S n (z) 

S 3 n (z) - zS 2 n (z) + (1 + r n {z))S n {z) - (1 + r n {z)f- = 0, z G C + , (11.8) 

where r n (z) := J R z{^pRz)-u From (11.1) and (6.4) we deduce, using Lemma 3.4, that 
(6.15) holds. Therefore, in view of (11.4), we obtain 

\r n {z)\ < ^p- < 53^= + — V q i(n)L qin < 5Arj ql (n)L qin < i- (11.9) 
Jn \/n n 10 



for z G C + . Note that we obtain the functional equation (11.8) from (6.17) replacing 

£ni(z) by r n (z) and e n2 (z) by -(1 + r n (z))z/n. 

11.3. The roots of the functional equation for S n (z). For every fixed z G C + consider 
the cubic equation 

P(z,w) := w 3 - zw 2 + (l + r n (z))w - (l + r n (z))- = 0. 

n 

As in Section 6 denote the roots of this equation by Wj = Wj(z), j = 1,2,3. Repeating 
the arguments of Subsection 6.4 we prove that 

z 10 3 
w 1 = - + f n (z), where |f n (z)| < — , (11.10) 
n n z 

and \vjj — z/n\ > 10 3 /n 2 , j = 2, 3, for z G D\. Hence w± ^ Wj for j = 2, 3 and z G D\. 

As in Subsection 6.5 we obtain that W2 ^ W3 for z G D4 := {z G C : < Qz < 3, \Rz\ < 

2 — hs}, where hs := C3 1 ^ e (r] q i(n)L qin + 1/n). Hence the roots W\(z),W2(z) and 1^3(2) are 

distinct for z G -D4. Moreover iui(z) satisfies, by (11.1) and (11.10), the inequality 

\wi(z)\<6/n, zeD 4 . (11-11) 

Using the arguments of Subsection 6.5 we deduce the formula (6.28) for z G D A where 
g(z) := (z—wi) 2 —A—Ar n (z)+Awi(z—Wi) 7^ 0, z G -D4. Then we rewrite the formula (6.28) 
as follows 

■=\{z + {-iy- W* 2 - 4 + f n (z)) - \ Wl (z), j = l,2, (11.12) 

where r n (z) := 2zw 1 (z) — 3w 2 (z) — Ar n (z). By (11.9) and (11.11), this function admits 
the bound, for z G D4, 

\f n (z)\ < 10\ Wl (z)\ + 3\ Wl (z)\ 2 + 4\r n (z)\ < 280 (p ql (n)L qin + 1/n) . (11.13) 

In the same way as in Subsection 6.5 we obtain that S n (z) = w 3 (z), z G -D4. Denote 
B 4 :=[-2 + /i 3 ,2-/i 3 ]. 
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11.4. Estimate of the integral J B \G^ w (x + ie) — G^ n (x + ie)\ dx for < e < 1. We 
obtain an estimate of this integral, using the inequality 

/ \G» w (x + ie) -G^(x + ie)\dx < / \G^ w (x + ie) - G/x n {x + ie)\ dx+ 

I \Gjx n (x + ie) - G^ n (x + ie)\dx. (11.14) 

J B A 

Evaluating the function G^ w {z) — G^ n (z) for z G D 4 in the same way as in Subsection 6.6, 
we arrive at the bound 

/ \G^ w (x + ie) - GfX n (x + ie)\dx < c(r] ql (n)L qin + l/n). (11.15) 
Now we conclude from (11.7) that 

G ^)-G^)= 1 -^^> *GC+. (11.16) 
Since | *SVi(^) j > 1/3 for z G C + , we see from (11.9) and (11.16) that 

JjG lt .(x + ie)-G fin (x + ie)\dx< J ^g^Vie)? ** ~ V' £G( °' 1] - (1L17) 
From (11.14), (11.15) and (11.17) we finally obtain 

J \G^ w (x + ie)-G^x + ie)\dx<c(^ ql (n)L qin + l/nj, e G (0, 1]. (11.18) 

11.5. Completion of the proof of Theorem 2.1. Note that 

/ p^(x)dx>l-hl /2 . (11.19) 
From (11.18) and (11.19) we deduce, using the Stieltjes-Perron inversion formula, 

<(£ 4 ) > l-c(^i(n)L, in + l/n). (11.20) 
Finally we deduce from (11.18)- (11.20) and the Stieltjes-Perron inversion formula that 

A(/<,/v) <c{rj ql (n)L qin + l/nj. (11.21) 
The statement of the theorem follows immediately from (11.5), (11.6) and (11.21). □ 

Proof of Corollary 2.2. The inequality (2.18) follows immediately from (2.17) and 
the Lyapunov inequality 1 = m^ 2 < fi]J q for q > 2. □ 
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